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PREFACE 


This  final  report  summarizes  a  fourteen  month  program  of  infrared 
data  collection  and  infrared  multispectral  target  cueing  analysis.  The 
research  was  carried  out  for  the  Night  Vision  Laboratory  (NVl)  of  the 
Army’s  Mobility  Equipment  Research  and  Development  Command  (MERDC)  at 
Ft.  Belvoir.  The  basic  objectives  of  this  effort  were  to  generate  a 
multispectral  infrared  data  base  of  tactical  targets  against  real  back¬ 
grounds  in  the  1-14  pm  region  and  to  use  these  data  to  develop,  analyze, 
and  evaluate  multispectral  target  cueing  algorithms.  The  motivation  behind 
the  effort  lay  in  the  potential  aid  to  new  generation  imagers  that  infrared 
multispectral  pattern  recognition  technology  may  add  to  already  developed 
spatial  and  3-D  cueing  techniques.  The  scope  of  the  program  included: 
hardware  modification  of  an  optical/mechanical  line  scanner  to  provide  high 
spatial  resolution  imagery;  test  plan  development  and  scanner  performance 
characterization,  airborne  infrared  multispectral  data  acquisition  at  Fort 
A.P.  Hill;  digital  reformatting  of  the  raw  data  to  computer  compatible 
tapes;  and  initiation  of  limited  analysis  of  the  data  to  investigate 
the  utility  of  multispectral  classification  cueing  techniques. 

The  research  described  in  this  report  was  performed  under  Contract 
DAAK70-77-C-0203  and  covers  the  period  22  September  1978  through  6  December 
1978.  Mr.  G.  David  Singer  served  as  Technical  Monitor  for  NVL.  The 
program  was  directed  by  Richard  R.  Legault,  Vice  President  of  the  Environ¬ 
mental  Research  Institute  of  Michigan  (ERIM)  and  Max  Bair,  Head  of  the 
Infrared  and  Optics  Division  Program  Office.  The  principal  investigator 
for  this  program  was  Stephen  R.  Stewart.  The  ERIM  number  for  this  report 
is  131500- 9-F. 

The  authors  wish  to  acknowledge  the  technical  contributions  provided 
by  J.  Ladd,  W.  Juodawlkis,  and  H.  Doss  in  the  data  collection  effort  and 
K.  Nalepka  and  R.  Horvath  for  their  valuable- comments  and  suggestions  during 
the  writing  of  this  report.  The  authors  also  thank  C.  Conley  and  F.  Elliott 
for  their  help  in  preparation  of  reports. 
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INTRODUCTION 


Studies  conducted  over  several  years  by  the  Environmental  Research 
Institute  of  Michigan  (ERIM)  and  others  have  contributed  to  the  develop¬ 
ment  of  advanced  spatial  and  3-dimensional  target  cueing  techniques. 

These  improved  image  processing  methods  have  reached  the  stage  where 
inclusion  of  automatic  cueing  technology  in  new  generation  imagers  now 
appears  feasible.  The  Army's  Night  Vision  Laboratory  (NVL)  at  Ft.  Belvoir 
has  been  assessing  the  state-of-the-art  of  these  currerit  techniques  through 
several  image  processing  study  contracts.  These  studies  have  indicated 
that,  although  thermal  (8-14  pm)  spatial  cueing  algorithms  and  multi- 
spectral  (MS)  pattern  recognition  techniques  in  the  visible  and  near- inf rared 
region  (.4  to  2.6  pm)  have  been  developed  which  Independently  provide 
successful  target  cueing,  none  have  considered  the  utility  of  the  entire 
infrared  spectrum  from  1  to  14  pm  simultaneously.  Not  only  was  there  a 
void  in  the  target  cueing  technology  development  in  this  respect,  a 
suitable  data  base  sufficient  for  a  thorough  study  did  not  exist. 

Because  of  the  potential  of  multispectral  pattern  recognition  as 
applied  to  the  entire  1  to  14  pm  spectral  region  to  supplement  other 
cueing  techniques  already  under  investigation,  NVL  and  ERIM  undertook 
a  two-phase  program  which: 

•  Provided,  in  the  1  to  14  pm  region,  an  airborne 
MS  data  base  of  tactical  targets  against  real 
backgrounds  from  which  MS  cueing  techniques 
could  be  evaluated;  and 

•  Performed  a  cursory  analysis  of  the  above 
data  to  determine  the  potential  utility  of 
classical  MS  techniques  for  thermal  target 
cueing. 

This  final  report  presents  the  results  of  that  effort. 

Phase  I  consisted  of  test  plan  development,  scanner  hardware  modifi¬ 
cations,  data  collection,  and  data  reformatting  and  comprised  most  of  the 
total  program  effort.  Since  this  infrared  data  base  was  intended  to  be 
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used  in  the  evaluation  of  various  MS  cueing  techniques,  an  emphasis  was 
placed  on  the  acquisition  of  a  complete  and  representative  data  matrix 
(in  terms  of  time  of  data  collection,  classes  of  tactical  targets,  and 
backgrounds/clutter  variability)  as  well  as  a  comprehensive  characteriza¬ 
tion  of  the  scanner  used  to  collect  the  data.  Section  2  of  this  report 
presents  a  brief  description  of  the  infrared  data  base  collected,  including 
an  outline  of  the  mission  plan  and  data  collection  scenario  as  well  as 
documentation  of  the  data  matrix.  More  detailed  information  concerning 
the  data  collection  effort  has  been  reported  under  separate  cover  [1], 

The  performance  of  ERIM's  M-7  multispectral  scanner  is  documented  in 
Section  3.  Only  those  performance  figures  which  are  unique  to  this 
mission  and  pertinent  to  the  analysis  of  the  imagery  are  reported.  A 
general  description  of  the  scanner's  operation  is  given  in  Appendix  A. 
Additional  performance  tests  and  a  discussion  of  techniques  for  radiometric 
calibration  of  the  data  are  described  in  Appendices  b  and  C  respectively. 

Phase  II  of  the  program  consisted  of  analysis  of  the  MS  data.  The 
emphasis  placed  on  collection  of  the  data  base  itself  severely  limited 
processing  and  subsequent  evaluation  of  classical  multispectral  "training 
set"  classification  techniques.  The  results  of  the  analyses  accom¬ 
plished  are  presented  in  Section  4.  There  are,  however,  several  other 
techniques  (such  as  multispectral  clustering,  pulse  height  and  pulse 
width  analysis,  and  MS  classification  combined  with  spatial  filtering) 
which  offer  potential  utility  in  target  cueing  and  which  should  be 
evaluated  against  these  data.  Recommendation  for  additional  target  cueing 
technique  evaluation  as  well  as  a  summary  of  the  conclusions  drawn 
from  the  limited  MS  analysis  is  reported  in  Section  5. 


[1]  Stewart,  Steve,  MS  Cueing  Techniques  Program,  Mission  Report  for 
Night  Vision  Laboratory  Program,  Camp  A.P.  "ill  Virginia  Test 
Site,  March  28,  29,  30,  1978,  Environmental  Research  Institute  of 
Michigan,  Ann  Arbor,  Mi.,  Report  No.  131500-8-P,  May  1978. 
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2.0 


DESCRIPTION  OF  INFRARED  MS  DATA  BASE 

2.1  INTRODUCTION 

Multlspectral  data  obtained  to  date  have  been  used  primarily  to 
evaluate  multlspectral  pattern  recognition  technology  applied  to  agri¬ 
cultural,  mineral  exploration,  water  quality  and  land  use  problems.  For 
such  applications,  night-time  operation  is  typically  not  required,  and 
targets  of  interest  to  tactical  military  scenarios  are  not  involved. 
Therefore,  virtually  all  multlspectral  data  collection  efforts  thus  far 
have  been  directed  at  the  evaluation  of  multlspectral  pattern  recognition 
techniques  in  non-military  applications  using  the  reflective  (0.3  to  3  pm) 
portion  of  the  spectrum.  Hence,  there  has  been  a  basic  deficiency  of 
multlspectral  imagery  involving  military  targets  in  tactical  military 
environments. 

The  deficiency  goes  farther  than  just  the  lack  of  targets  and  back¬ 
grounds,  however.  For  military  applications,  there  are  potential  target 
cues  which  occur  in  the  mid  and  far  infrared  region  (3  to  14  pm)  which 
have  not  been  exploited  with  today's  available  multlspectral  imagers. 

At  wavelengths  beyond  3  pm,  the  primary  source  of  infrared  energy  is  the 
thermal  self-radiation  produced  by  the  target  and  it  is  the  fact  that 
this  IR  radiant  energy  source  is  internalized  to  the  target  that  opens  a 
new  dimension  of  possibilities  for  multlspectral  pattern  recognition 
techniques. 

For  example,  the  identification,  composition,  and  material  properties 
of  both  man-made  and  natural  scene  objects  can  often  be  inferred  by  applying 
qualitative  knowledge  of  how  different  objects  respond  thermally  to 
various  environmental  driving  forces.  Since  the  radiometric  temperature 
of  an  object  depends  upon  its  surface  optical  characteristics  and  tempera¬ 
ture  and  since  different  materials  possess  different  optical  characteristics 
and  thermal  properties,  material  type  and  structural  differences  are  often 
Indicated  by  different  temperature  responses  under  equivalent  environmental 
driving  forces.  In  addition  these  material  properties  (e.g.  emissivity) 
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are  a  function  of  wavelength  and  thus  objects  of  the  same  surface  tem¬ 
perature  but  with  a  spectrally  different  material  properties  can  exhibit 
completely  different  IR  radiometric  temperatures  when  observed  in  different 
spectral  regions.  It  is  this  variation  in  radiometric  response  as  a 
function  of  material  type,  spectral  IR  properties,  and  a  given  environ¬ 
mental  condition  that  can  be  exploited  by  multispectral  pattern  recognition 
techniques  in  the  infrared  region. 

There  are  other  IR  cues  which  have  yet  to  be  fully  exploited  in 
the  military  application.  Consider  the  fact  that  the  temperature  of  an 
object  responding  only  to  normal  environmental  driving  forces  will  be 
different  from  that  of  a  like  object  within  which  heat  is  being  generated, 
or  upon  which  anomalous  heat  is  externally  applied.  Thus,  a  motor  casing 
is  hotter  when  the  motor  runs  (internal  heat  generation),  and  an  exhaust 
deflector  is  hotter  when  exhaust  is  present  (external  heat  application). 

It  is  a  fact  that  some  targets  are  recognizable  by  virtue  only  of  their 
"hot  spot"  distribution  while  others  have  "warmer  areas"  on  the  target 
surface  rather  than  hot  spots  because  of  conductive  or  convective  heat 
redistribution.  Infrared  image  exploitation  of  this  nature  utilizes  the 
qualitative  spatial  characteristics  of  objects  within  a  scene.  Only 
the  existence  of  contrasts  and  not  their  temperature  inferences  is  directly 
relevant  to  this  spatial  exploitation.  Hence,  night-time  IR  imagery, 
with  no  requirement  for  natural  illumination  and  with  familiar  spatial 
representation  of  targets  can  substitute  for  visible  Imagery  under  low 
light  level  conditions. 

Also,  infrared  shadows,  while  not  enjoying  a  direct  spatial  relation¬ 
ship  with  a  target,  can  still  provide  useful  cueing  information.  The 
magnitude  of  the  IR  shadow  contrast  is  a  function  of  how  long  the  target 
or  background  object  producing  the  shadow  is  stationary.  This  is  one 
case  wherein  both  the  thermal  IR  and  reflective  signals  are  positively 
correlated.  When  used  in  conjunction  with  spatial  processing  the  length 
and  width  of  the  shadow  may  provide  information  about  target  h  ight  or 
provide  other  inferences  about  the  object  producing  the  shadow. 
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Investigation  of  the  utility  of  the  IR  target  cues  suggested  in  the 
previous  paragraphs  in  discriminating  military  vehicular  targets  from 
background  in  a  tactical  military  environment  requires  a  data  base  of 
multispectral  imagery  of  the  appropriate  targets  and  background  which 
can  be  used  to  evaluate  not  only  the  significance  of  the  spectral  signa¬ 
tures  but  also  the  variation  in  scene  spatial  distribution  encountered  under 
a  variety  of  environmental  conditions,  vehicle  types  and  operating  modes, 
and  background  types.  The  generation  of  a  suitable  data  matrix  which 
takes  into  consideration  all  these  factors  was  considered  a  major 
responsibility  of  the  data  collection  effort.  The  salient  features  of 
the  data  matrix  are  discussed  in  the  following  sections. 

2.2  SPECTRAL  CONSIDERATIONS 

Paramount  to  the  successful  application  of  any  multispectral  cueing 
technique  (visible  or  IR)  is  the  choice  of  wavelengths  to  use  and  the 
associated  spectral  differences  in  target /background  optical  and  thermal 
characteristics  encountered  in  the  chosen  wavelength  regions.  Obviously, 
if  wavelength  regions  are  chosen  where  little  difference  in  target  and  back¬ 
ground  optical  cueing  properties  exist,  then  any  multispectral  cueing  tech¬ 
nique  based  on  spectral  signatures  alone,  is  doomed  to  failure.  In 
like  manner,  if  the  regions  are  chosen  without  regard  for  1)  the  effect 
of  spectral  absorption  and  scattering  within  the  intervening  atmosphere 
and  2)  the  detector  spectral  operating  characteristics  then  differences 
in  target/background  optical  properties  may  be  masked  by  noise.  Hence, 
the  spectral  band  requirements  for  the  multispectral  data  base  were 
based  upon  the  following  criteria: 

(1)  Atmospheric  Transmission 

(2)  Target /Background  Optical  Properties 

(3)  Detector  Noise  Considerations 

As  shown  in  Figure  1,  there  are  several  broadband  regions  in  the  1 
to  14  pm  region  that  can  be  considered  from  an  atmospheric  transmission 
standpoint.  These  are: 
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TRANSMITTANCE  (%) 


ATMOSPHERIC  TRANSMISSION 


Transmission  of  1000  foot  horizontal  air  path  at  sea  level, 
5.7  mm  precipitable  wafer,  79T  temperature. 
(Reprinted  with  permission  of  Santa  Barbara  Research  Center) 
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(a) 

1 

to 

1.1  pm 

(b) 

1. 

,1 

to  1.4 

pm 

(c) 

1. 

,5 

to  1.8 

pm 

(d) 

2 

to 

2.6  pm 

(e) 

3 

to 

4.2  pm 

(f) 

4. 

.5 

to  5.5 

Mm 

(g) 

8 

to 

14  pm 

spectral  bands. 

Any  application  (military  or  other)  of  multispectral  processing  techniques 
to  infrared  imagery  of  the  earth  viewed  through  the  atmosphere  in  the 
1  to  14  Pm  spectral  region  is  limited  to  these  broad  spectral  bands; 
outside  of  these  bands,  the  atmosphere  becomes  the  prime  source  of  the 
IR  radiant  energy  observed  in  the  imagery. 

Of  these  seven  (7)  broadbands,  the  first  four  are  considered  reflective 
bands  and  the  last  two  emissive  bands.  Figure  2  shows  the  spectral  depen¬ 
dence  of  the  reflected  solar  and  self-emitted  radiant  energy  for  an  object 
with  a  spectrally  constant  reflectance  of  20%,  an  emissivity  of  80%  and 
no  atmospheric  effects.  As  can  be  seen,  the  spectral  region  where  these 
two  radiant  energy  components  are  equal  is  in  the  3  to  4  pm  region.  Hence, 
band  five  as  defined  by  the  atmospheric  limits  above  can  be  classified 
neither  totally  reflective  nor  totally  emissive  as  though  in  the  absence 
of  any  solar  energy.  The  target  radiance  observed  in  this  band  (though 
admittedly  low  in  magnitude)  is  primarily  self -emission.  Figure  3  shows 
this  effect  more  clearly.  Spectral  apparent  temperature  measurements  in 
the  3.3  to  14  pm  region  were  obtained  from  bare  soil  on  a  clear  (.1 
cloud  cover)  and  and  overcast  day  [Ref.  2].  The  increase  in  apparent 
soil  temperature  in  the  3  to  4  pm  region  on  the  clear  day  (top  graph) 
is  due  entirely  to  reflected  solar  energy.  No  increase  is  observed  in  this 


[2]  Spectral  and  Polarization  Characteristics  of  Selected  Targets  and 
Backgrounds:  Instrumentation  and  Measured  Results  (3.3  -  14.0  pm), 
D.  Faulkner,  R.  Horvath,  J.P.  Ulrich,  E.  Work,  Willow  Run  Labora¬ 
tories,  The  University  of  Michigan,  Technical  Report  No.  AFAL-TR- 
71-199,  August  1971. 
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Sunlight  Irradiance,  Thermal  Emission  (W/cm  um) 


! 


(Wavelength  (urn) 

FIGURE  2.  SOLAR  ( 6000*  K)  REFLECTED  AND  THERMALLY  EMITTED 
RADIANCE  ( ^=0.8,  *  =0.2)  NEGLECTING  ATM  OSPHERIC 
EFFECTS 


18 


FIGURE  3.  BARE  SOIL  AND  SKY  SPECTRAL  APPARENT  TEMPERATURE 
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region  on  the  overcast  day  since  the  sun  is  obscured  by  clouds  and  the 
only  energy  reflected  from  these  in  the  3  to  4  urn  region  is  that  coming 
from  the  cloud  layer. 

It  is  obvious  from  this  discussion  that  a  single  detector  band  covering 
the  3  to  4  pm  region  could  provide  ambiguous  thermal  target  cues  when 
operated  alone  in  clear  day  conditions.  An  anamolously  high  apparent 
scene  object  temperature  obtained  with  such  a  detector  could  be  due  to 
either  a  hot  (running)  vehicle  or  a  solar  glint  from  some  naturally 
occuring  smooth  surface  like  water.  If  this  band  is  to  be  used  in  day¬ 
time  scenarios,  then  this  ambiguity  must  be  removed.  The  most  obvious 
way  is  to  check  for  the  presence  of  the  solar  glint  in  the  signal  by 
looking  for  it  in  a  band  that  is  totally  reflective.  Thus,  if  the  same 
high  target  signal  as  observed  in  the  3  to  4.2  pm  band  is  simultaneously 
observed  in  a  band  somewhere  in  the  0.4  to  3  pm  region,  then  a  solar 
glint  is  being  observed  and  not  necessarily  a  'hot''  target. 

Given  the  broadband  limits  defined  previously,  it  is  likely  that 
some  spectral  differences  in  target  and  background  optical  characteristics 
and  thermal  properties  exist  not  only  across  the  listed  broadbands  but 
also  within  a  given  band.  The  questions  remain  however:  1)  whether 
the  observed  radiometric  temperature  difference  associated  with  the 
material  property  differences  are  large  enough  to  observe  with  a  multi- 
spectral  scanner  and  2)  whether  the  measured  radiometric  difference  can 
be  positively  correlated  to  a  specific  target  and  background. 

Figures  4  and  5  present  emittance  (e)  and  reflectance  (Ssl  -  e) 
curves  in  the  1  to  15  pm  region  for  several  backgrounds  and  one  O.D. 
military  paint.  As  can  be  seen,  most  of  the  curves  within  the  defined 
atmospheric  windows  are  smoothly  varying  functions.  The  one  obvious 
exception  occurs  at  between  9  and  10  pm  in  the  O.D.  paint  data  (Figure  4) 
where  a  narrow  reststrahlen  band  is  observed  due  to  the  presence  of  SK>2 
in  the  paint.  Unfortunately,  this  Si02  reststrahlen  band  is  also  present 
in  soil  (look  around  9  pm  in  Figure  3)  and  as  Figure  3  indicates,  the 
magnitude  is  dependent  upon  the  degree  of  cloud  cover.  Because  of  this 
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Pressed  Dormant  Maple  Leaf  (Dry,  Top) 


Sand  No.  113  (Maury  Silt  Loam,  Tennessee) 


Sand  No.  115  (Dublin  Clay  Loam,  California) 


FIGURE  5.  SPECTRAL  REFLECTANCE  OF  SIX  TERRAIN 
FEATURES  IN  THE  1-15  pm  REGION. 
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environmental  dependence  and  because  many  military  paints  do  not  contain 
SK^,  there  is  no  unambiguous  technique  employing  a  SiOj  reststrahlen  band 
detector  to  determine  the  presence  of  painted  metal.  Other  than  the 
reststrahlen  band  phenomena,  military  paints  in  general  do  not  show 
rapidly  varying  spectral  reflectance  and/or  spectral  emittance  charac¬ 
teristics  in  the  1  to  14  pm  spectral  region  [3],  Because  of  the  above 
factors  and  adding  the  additional  fact  that,  except  for  the  8  to  14  pm 
spectral  band,  consideration  of  more  than  two  independent  sub-bands  within 
each  of  the  atmospheric  window  regions  is  precluded  by  detector  signal- 
to-noise  considerations,  multiple  sub-bands  within  any  of  the  broad 
atmospheric  windows  were  not  considered  for  the  proposed  effort. 

Shown  in  Table  1  are  the  spectral  bands  chosen  for  this  MS  data  base 
along  with  a  summary  of  each  band's  utility.  The  first  band  in  Table  1 
corresponds  to  the  1  to  1.1  pm  band  described  earlier.  Detector  signal 
to  noise  considerations  dictated  a  0.8  pm  cut-on  instead  of  1  pm.  The 
1.1  to  1.4  pm  and  1.5  to  1.8  pm  bands  were  not  included  because  of  detector 
hardware  limitations  in  the  collecting  sensor.  However,  significant 
amounts  of  background  data  have  been  collected  in  these  bands  in  recent 
years  on  non-military  programs  although  no  known  vehicular  data  exists. 

2.3  TARGET  AND  BACKGROUND  CONSIDERATION 

The  choice  of  spectral  bands  for  the  multispectral  data  collection 
and  cueing  evaluation  effort  was  made  somewhat  objectively.  However, 
consideration  of  targets  and  backgrounds  was  much  more  subjective  in  that 
many  different  tactical  environments  can  be  presented  to  a  sensor,  depending 
upon  battlefield  scenario  and  season.  Obviously,  not  all  target  vehicle 
and  background  combinations  need  be  considered  since  some  combinations 
occur  too  infrequently  to  warrant  significant  effort  or  the  combination 
lends  itself  to  conventional,  visible  cueing  techniques  (e.g.  O.D.  painted 


[3]  Handbook  of  Military  Infrared  Technology,  William  L.  Wolfe,  Editor, 
The  University  of  Michigan,  1965,  Sponsored  by  Office  of  Naval 
Research,  Department  of  the  Navy,  Washington,  D.C. 


23 


TABLE  1.  NOMINAL  MS  DATA  BASE  SPECTRAL  BANDS 


Channel 


Nominal  Band  Limits 


Solar  glint  discriminator 
plus  reflective  multispectral 


Solar  glint  discriminator 
plus  reflective  multispectral 


Reflective  multispectral, 
emissive  multispectral,  and 
thermal  temperature 


Thermal  temperature  plus 
emissive  multispectral 


tank  in  the  desert).  Even  without  these,  however,  the  number  of  combina¬ 
tions  remaining  is  still  much  too  large  to  attempt  a  data  collection 
effort  involving  all  possibilities.  Hence,  for  this  particular  program 
a  decision  was  made  in  the  interest  of  time  and  money,  to  limit  the 
target/background  combinations  to  one  which  was  of  the  most  interest  of 
the  sponsor  —  namely,  an  Eastern  European-like  scenario. 

The  site  selected  by  the  sponsor  for  this  experiment  was  the  test 
facility  located  at  Camp  A.P.  Hill,  Virginia.  The  backgrounds  available 
at  the  time  of  data  collection  flights  were  primarily  baresoil  in  open 
fields,  tree  stands  of  leafless  deciduous  and  coniferous  trees  (types 
unknown)  and  forest  carpet  of  leaves,  broken  limbs,  and  soil.  Although 
this  site's  latitude  is  approximately  15°  south  of  the  latitude  of  Germany, 
variety  of  backgrounds  encountered  here  were  thought  to  be  reasonably 
representative  of  what  might  be  found  in  some  spring-time  European 
scenarios . 

This  was  based  in  part  on  the  presumption  that  pine  tree  canopies, 
leafless  tree  trunks,  and  spring-time  forest  carpets  presumably  don't 
exhibit  significantly  different  thermal  spectrums  over  the  many  varieties 
in  existence  at  these  latitudes.  While  it  is  quite  likely  that  more 
northernly  latitudes  do  have  different  soil  mixtures  and  types  of  hard¬ 
wood  and  coniferous  tree  stands  than  present  at  this  site,  it  is  equally 
likely  that  variations  exist  at  a  given  latitude  as  a  function  of  longitude 
so  that  one  may  question  how  good  a  European  simulation  any  U.S.  site 
is.  Since  this  particular  site  had  the  "kinds''  of  spring-time  backgrounds 
expected  in  Europe  and  since  the  facilities  existed  to  supply  the  tactical 
targets  and  other  logistical  ground  support,  the  site  was  used.  (Figure 
6  shows  an  airborne  photograph  of  the  target  site.) 

Eleven  targets  were  made  available  for  the  data  collection  exercise. 
These  included  the  following: 

(1)  M-60  Tank  (3,  2  running  [1  diesel,  1  gas]  and 
1  static) 

(2)  M-113  APC  (2  diesel,  2  gas,  all  running) 


FIGURE  6.  AERIAL  PHOTOGRAPH  OF  TARGET  SITE(TARGETS  IN  POSITION) 

3/30/78,  PASS#5.  1830 


(3)  3/4,  2  1/2  and  5  Ton  Trucks  (all  running) 

(4)  Jeep  (running) 

The  targets  were  started  up  about  one-half  hour  before  the  overflight 
and  kept  running  throughout  the  period  of  data  collection.  The  targets 
were  positioned  in  one  of  the  twelve  target  locations  shown  in  Figure  7 
and,  after  each  data  collection  pass,  were  moved  to  the  next  available 
target  position.  This  permitted  each  class  of  vehicle  to  be  placed  in 
a  background  of  varying  clutter  and/or  cover  while  maintaining  a 
vehicle  operating  temperature  consistent  with  those  that  might  be  en¬ 
countered  in  an  operational  tactical  scenario.  A  summary  of  the  backgrounds 
at  each  target  location  is  shown  in  Table  2.  The  data  matrix  of  target 
positions  for  each  flight  is  shown  in  Tables  3  through  6. 

The  data  collection  mission  took  place  over  the  period  from  28  March 
through  1  April  1978.  ERIM's  M-7  scanner  and  associated  boresighted  camera 
were  the  principal  sensors  flown  in  our  C-47  aircraft.  The  instrumentation 
configuration  is  depicted  in  Figure  3.  Four  flights  were  planned  so  that 
data  could  be  collected  under  conditions  producing  maximum  and  minimum 
scene  contrasts:  0930,  1330,  1830,  and  2330  hours.  The  general  profile 
for  each  of  the  four  missions  consisted  of  a  total  of  ten  (1)  passes 
ranging  in  altitude  from  6400  ft.  to  800  ft.  This  provided  ground 
resolutions  for  the  normal  nadir  viewing  which  range  from  1.6  ft.  square 
to  12.8  ft.  square.  The  only  exception  to  this  occurred  during  the  0930 
flight  when  four  additional  passes  at  low  altitude  were  completed  which 
utilized  slightly  different  spectral  bands  and  a  forward  tilt  of  the 
scanner  of  55°  (35°  depression  angle).  These  passes  were  flown  to  provide 
data  which  could  be  used  for  additional  spectral  analysis  and  which  would 
represent  a  scene  aspect  more  typical  of  those  viewed  by  forward  looking 
air-to-ground  sensors.  Table  7  and  8  show  the  altitude,  time  of  day  and 
spectral  bands  for  each  of  the  passes  flown  during  the  four  flights. 

The  data  base  just  described,  consisting  of  data  from  all  forty-four 
(44)  data  collection  passes  over  the  A.P.  Hill  site,  is  now  resident 
on  a  set  of  14  computer  compatible  tapes  (CCTs).  Each  file  of  these 
tapes  contains  five  channels  of  spectral  data  as  described  in  Table  9. 
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TABLE  2.  DESCRIPTION  OF  BACKGROUND  AT  TARGET  LOCATIONS 


Target  Position  Background  Description 

0  Open  soil  mixed  with  sparsely  populated 

pine  trees.  Completely  open  on  one 
side.  Approximately  15%  obscuration. 

1  Open  soil  cut  out  of  dense  stand  of 

pine  trees.  Open  soil  on  one  side  — 
no  overhead  obscurations. 

2  Forest  carpet  surrounded  closely  on 

three  sides  by  dense  pine  tree  stand 

—  50%  overhead  obscuration. 

3  Forest  carpet  loosely  bounded  by 

dense  pine  tree  stand  on  one  side  and 
mixture  of  pine  trees  and  leafless 
deciduous  trees  on  other  —  no  over¬ 
head  obscurations. 

4  Narrow  dird  road  boardered  closely 

on  both  sides  by  dense  pine  trees  — 

5%  obscuration. 

5  Forest  carpet  loosely  surrounded  by 

leafless,  deciduous  trees  and  pine 
trees  —  no  overhead  obscurations. 

6  Forest  carpet  within  tree  stand  of  pines 

and  leafless  deciduous  -  15%  overhead 
obscuration. 

7  Forest  carpet  in  small  pine  tree  stand 

—  not  very  dense  —  25%  obscuration. 

8  Forest  carpet  in  sparsely  populated, 

leafless  deciduous  tree  stand  —  10% 
overhead  obscuration. 

9  Open  forest  carpet  cut  out  of  dense 

pine  tree  stand  —  no  obscuration 

Forest  carpet  loosely  located  between 
small  pine  tree  stand  on  one  side  and 
sparsely  populated  leafless  deciduous 
trees  on  other  —  no  obscuration. 
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TABLE  2.  DESCRIPTION  OF  BACKGROUND  AT  TARGET  LOCATIONS  (Continued) 


Target  Position 


Background  Description 

Open  soil  bounded  lossely  on  one 
side  by  small  stand  of  pine  trees 
—  no  obscuration. 


Open  soil  bounded  closely  on  one  side 
by  small  stand  of  pine  trees  — 


TABLE  3.  MATRIX  OF  VEHICLE  POSITIONS  DURING  3/28/78,  2330  HOURS 
DC-3  DATA  PASSES  OV'JR  SITE 


Cold  Tank 


TABLE  4.  MATRIX  OF  VEHICLE  POSITIONS  DURING  3/29/78,  1830  HOURS  FLIGHT 
DC-3  DATA  PASSES  OVER  SITE 
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TABLE  5.  MATRIX  OF  VEHICLE  POSITIONS  DURING  3/30/78,  0930  HOURS  FLIGHT 
DC-3  DATA  PASSES  OVER  SITE 


Cold  Tank 


TABLE  6.  MATRIX  OF  VEHICLE  POSITIONS  DURING  3/30/78,  1330  HOURS 
DC-3  DATA  PASSES  OVER  SITE 
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Note:  This  flight  was  started  with  the  vehicles  in  the  same  positions  as  the  final 
from  Flight  //3.  The  jeep  in  Position  If 2  began  to  leak  oil  on  the  first  run  / 
was  shut  off  and  moved  slightly  out  of  position  where  it  remained  for  the  re 
of  the  passes. 


TABLE 

7 .  ALTITUDES 

,  TIMES,  AND 

BANDS  DURING 

EVENING  MISSIONS 

Pass  # 

Alt.  (ft.) 

3/29/78 
Local  Time 

3/29/78 
Local  Time 

Spectral  Bands 

1 

6400 

0022 

1908 

4. 3-5. 3  pm  8.4-13.7  pm 

2 

3200 

0034 

1916 

••  fl 

3 

1600 

0041 

1922 

tl  •• 

4 

800 

0046 

1927 

•1  It 

5 

800 

0051 

1932 

••  •• 

6 

800 

0055 

1936 

••  II 

7 

800 

0059 

1947 

••  If 

8 

3200 

0107 

1950 

II  II 

9 

1600 

0114 

1956 

•1  II 

10 

800 

0119 

2002 

It  II 
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TABLE  8.  ALTITUDES,  TIMES,  AND  SPECTRAL  BANDS  DURING  DAYTIME  MISSIONS 


Pass  # 

Alt.  (ft.) 

3/30/78 
Local  Time 

1 

6400 

1025 

2 

3200 

1033 

3 

1600 

1039 

4 

800 

1043 

5 

800 

1047 

6 

800 

1051 

7 

800 

1055 

8 

3200 

1102 

9 

1600 

1108 

10 

800 

1112 

11* 

800 

1119 

12** 

500 

1137 

13** 

900 

1142 

14** 

900 

1149 

15** 

900 

1156 

3/29/78 

Spectral  Bands 

Local  Time 

(pm) 

1357 

.86-1.04;  1.9- 

2.55; 

2. 9-3. 9;  4.3 

-5.3; 

1403 

8.4-13. 

7; 

1408 

It  It  II 

•1 

1412 

II  II  II 

•I 

1418 

II  II  II 

•1 

1421 

•1  II  II 

II 

1425 

II  II  II 

II 

1431 

II  II  II 

•  1 

1436 

II  II  II 

II 

1440 

•1  II  II 

II 

II  II  II 

II 

II  II  II 

II 

II  II  II 

II 

.86-1.04;  3.5- 

4.0; 

3. 9-4. 7;  8.4 

-13.7 

.86-1.04;  2.0- 

2.6; 

3. 5-4.0;  3.9 

-4.7; 

*  Did  not  cover  target  site 

**  Scanner  tilted  35°  down  fron  horizon 


TABLE  9.  CCT  TAPE  CHANNEL- SPECTRAL  REGION  DOCUMENTATION 


Nominal  Spectral  Region  for  Each  Tape  Channel  (pm) 


Tape  No. 

Channel  1 

Channel  2 

Channel  3 

Channel  4 

Channel  5 

NVL001A 

.78-1.1 

2  -  2.6 

3  -  4.0 

4.5  -  5.5 

8  -  13.5 

2A 

•V 

•1 

II 

II 

II 

3A 

If 

If 

II 

II 

•1 

4A 

II 

ft 

II 

II 

•  I 

5 

ft 

If 

II 

II 

•1 

6 

If 

•1 

II 

II 

II 

7 

If 

II 

II 

II 

•I 

8 

If 

II 

II 

•1 

«« 

9 

II 

II 

It 

II 

II 

10 

If 

•1 

II 

II 

•1 

11  (Run  14)  " 

- 

3.9  -  4.7 

3.5  -  3.9 

II 

11  (Run  15)  " 

2  -  2.6 

ft 

II 

- 

12 

If 

ii 

3  -  4.0 

4.5  -  5.5 

8  -  13.5 

13 

II 

ii 

II 

II 

II 

14 

If 

ii 

It 

II 

•I 

The  CCTs  contain  132  images  of  M-60  tank-types,  176  Images  of  M-113 
APC-type,  176  Images  of  truck-type  targets.  The  flight  parameters 
associated  with  each  of  the  aircraft  overpasses  are  described  in 
detail  in  Reference  [4.4]. 
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Further  documentation  of  the  organization  of  the  data  base  is  provided 
in  Table  10.  For  those  interested  in  using  this  data  base  to  test  target 
cueing  techniques,  the  tape  format  is  described  in  Appendix  D. 
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TABLE  10.  DOCUMENTATION  FOR  COMPUTER  COMPATIBLE  TAPES 
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line  number  (maximum  or 


3.0 


M-7  SYSTEM  PERFORMANCE  CHARACTERIZATION 

In  order  that  the  information  contained  in  the  data  base  be  maxi¬ 
mally  useful,  it  is  important  that  the  user  be  familiar  with  the 
performance  characteristics  of  the  sensor  system.  We  provide  discussions 
of  these  characteristics  in  this  section.  This  information,  along  with 
that  contained  in  the  appendices,  completely  documents  the  operation  and 
performance  of  the  sensor  employed  to  generate  the  data  base. 

3.1  BACKGROUND 

The  airborne  line  scanner  used  on  the  data  collection  task  was  modi¬ 
fied  to  meet  the  requirements  specified  for  the  data  set  by  the  Sponsor. 
The  modifications  included  the  following: 

(1)  Increasing  the  scan  rate  of  the  scanner  to 
obtain  contiguous  ground  coverage  for 
minimum  ground  resolution  (.5  m) ; 

(2)  The  purchase  and  installation  of  a 
3-element  detector  array  to  obtain 
spatially  registered  multispectral  data 
for  the  near  IR  (2.0  to  2.6  pm)  and  mid 
IR  (3.0  to  5.5  pm)  spectral  regions;  and 

(3)  Interfacing  a  high  data  rate  digital  recorder 
to  the  modified  M-7  scanner. 

The  results  of  these  modifications  were  to  alter  the  M-7's  performance 
characteristics  from  its  nominal  configuration,  a  general  description 
of  which  is  given  in  Appendix  A. 

In  the  following  text,  the  results  of  the  laboratory  tests  documenting 
the  performance  of  the  M-7  scanner  configured  for  the  A.P.  Hill  mission 
are  summarized.  Only  a  limited  discussion  is  given  on  the  majority  of 
procedures  used.  Further  detail  on  the  M-7  scanner  and  some  of  the 
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characterization  tests  can  be  found  in  Appendix  B.  Nominal  operating 
characteristics  of  the  system  as  it  was  flown  for  the  program  are 
shown  in  Table  11. 

3.2  SYSTEM  SPECTRAL  RESP0NS1VITY 

Figures  8  through  11  show  the  actual  relative  spectral  response  of 
the  five  detector /filter  combinations  used  in  the  MVL  configuration. 
Spectral  data  for  the  near  infrared  band  was  obtained  by  actual  measure¬ 
ment  of  the  M-7  spectrometer  using  a  monochroneter  at  the  focal  plane  of 
a  laboratory  collimator.  The  mid-lR  spectral  responsibities  were  measured 
by  the  vendor  supplying  the  InSb  3-element  array.  A  curve  for  the  long 
wave  IR  band  was  generated  from  the  detector  spectral  responsibity  measured 
by  the  detector  vendor  and  a  spectral  cut-on  filter  transmission  measured 
at  ERIM.  Except  for  the  near  IR  band  which  was  measured  through  the  entire 
scanner,  these  data  represent  the  spectral  response  of  only  two  elements  in 
the  optical  chain,  i.e.,  narrow  band  spectral  filter  and  the  detector. 
However,  since  these  are  the  major  components  in  determining  spectral 
responslvity,  the  other  components  in  the  optical  path  such  as  the  scanning 
mirror  and  focusing  optics  can  be  neglected  without  introducing  significant 
error. 

The  50%  spectral  response  points  as  depicted  in  Figures  4  through  6 
are  listed  below: 


(1) 

.86  - 

1.04 

(2) 

1.9  - 

2.55 

(3) 

2.9  - 

3.9 

(4) 

4.3  - 

5.3 

(5) 

8.4  - 

13.7 

3.3  SPATIAL  REGISTRATION 

Figures  12  and  13  show  the  spatial  registration  of  the  five  (5)  spectral 
channels  by  viewing  a  point  target  of  appeoximately  2  1/2  mr  size  with  the 
M-7  scanner  placed  on  a  laboratory  collimator.  As  the  figures  show,  two 
of  the  channels  are  not  in  registration  in  the  cross -track  dimension.  These 
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TABLE  11.  M-7  SCANNER  PERFORMANCE  CHARACTERISTICS 

(A.P.  Hill  Configuration) 

•  S  Spectral  Bands  in  the  Near,  Mid,  and  Far  IR  Regions 

•  90°  External  FOV  (+  45°  from  nadir) 

•  2.0  mrad  Spatial  Resolution 

•  0.2°C  Nominal  (mid  and  far  IR)  Thermal  Resolution 

•  5  Inch  Diameter  Collector  Optics 

•  100  Scans /Sec 

•  DC  to  150  KHz  Electonlc  Bandwidth 

•  Roll-Stabilized  Imagery 

•  8-Bit  Airborne  High  Density  Digital  Recording  (10  K  bits/ inch) 
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Wavelength  (jim) 

FIGURE  9.  SPECTRAL  RESPONSE  OF  LONG  WAVE  IR  DETECTOR/ FILTER  COMBINATION  (HgCdTe) 


FIGURE  10.  SPECTRAL  RESPONSE  OF  MID-IR  InSb  DETECTOR  ARRAY 


FIGURE  12.  SPATIAL  REGISTRATION  OF  NEAR  -IR  BAND 


(.86-1.04  um)  WITH  THE  MID-IR  BANDS 
(1.9-2.55, 2.9-3.9,-nd  4.3-5.3  urn) 


FIGURE  13.  SPATIAL  REGISTRATION  OF  NEAR-IR  BAND 
(.86-1.04  fjm)  WITH  THE  THERMAL  BAND 
(8.4-13.7  /jm) 


are  two  channels  of  the  three-element  InSb  array  and,  because  of  the 
construction  of  the  array,  these  channels  are  r.  :t  expected  to  be  in 
registration  in  this  dimension.  However,  the  two  offset  channels  were 
brought  into  cross-track  registration  (within  1/2  pixel)  with  the  other 
three  channels  by  digital  deskewing  during  the  HDDT  to  CCT  reformating 
process.  The  registration  in  the  in-track  dimension  for  all  channels  is 
within  1/4  pixel.  All  data  produced  in  CCT  format  on  this  program  is 
spatially  registered  to  within  one-half  pixel. 

3.4  INSTANTANEOUS  FIELD-OF-VIEW 

The  instantaneous  f ield-of-view  or  resolution  of  the  M-7  system  as 
configured  for  this  data  collection  mission  was  calculated  to  be  nominally 
2.0  mr  for  all  channels.  However,  the  most  useful  measure  of  the  resolu¬ 
tion  of  an  imaging  system  is  the  MTF  (modulation  transfer  function)  which 
represents  the  sine-wave  spatial  frequency  amplitude  response  of  the 
system.  The  prime  utility  of  the  MTF  as  a  measure  is  that  it  permits 
the  cascading  of  effects  of  several  major  components  in  determining  the 
resolution  of  the  overall  system.  At  a  given  spatial  frequency,  the 
overall  MTF  is  simply  the  product  of  the  MTF's  of  all  components.  The 
MTF's  presented  in  Figures  14,  15,  16,  17  and  18  represent  the  overall 
system  MTF  for  each  channel.  The  spatial  frequency  is  given  in  cycles/ 
milliradians. 

Examination  of  these  figures  point  out  several  important  facts. 

First,  it  can  be  seen  that  in  general  the  cross-track  MTF  is  less  than 
* 

the  in-track  MTF  .  This  is  to  be  expected  to  some  degree  since  the 
response  of  signals  at  the  higher  spatial  frequencies  in  the  cross-track 
dimension  is  limited  somewhat  by  detector  time  constant  and  amplifiying 
electronics  bandwidth.  The  high  spatial  frequencies  in  the  in-track 
dimension,  being  sampled  at  a  much  lower  electronic  frequency  defined 

* 

By  cross-track  we  mean  that  dimension  which  is  perpendicular  to  the 
flight  line  direction  and  in  the  direction  of  mirror  rotation.  The  in¬ 
track  dimension  is  in  the  direction  of  the  aircraft  flight  line. 
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Detector:  PM  9-1 
Wavelength:  0.78  -  1.11 


FIGURE  14.  MTF  VS  SPATIAL  FREQUENCY 


FIGURE  15.  MTF  VS  SPATIAL  FPEQUENCY 


tector:  3-Element  InSb  3-8 


FIGURE  17.  MTF  VS  SPATIAL  FREQUENCY 


by  the  mirror  rotation  rate,  is  not  limited  by  these  same  factors.  The 
general  trend  of  Figures  9  through  13  is  one  in  which  the  difference 
between  cross-track  and  in-track  MTF  increases  as  the  detectors  go 
from  photomultiplier  tube  (Figure  14)  to  HgCdTe  (Figure  13)  to  Indium 
Antimonde  (Figures  10,  11,  and  12)  is  indicative  of  the  detector  response 
time  influence  on  cross-track  MTF.  A  photomultiplier  tube  has  the  fastest 
response  time  followed  by  the  HgCdTe  and  finally  the  InSb  detectors. 

The  "dip"  in  response  at  some  mid-frequencies  ('v-  0.15  cycles/mr) 
in  Figures  15,  16,  and  17,  however,  is  not  related  to  any  electronic 
frequency  response  problem.  Rather,  it  is  caused  by  undesirable  surface 
reflections  which  occur  because  of  manufacturing  defects  in  the  three- 
element  InSb  detector  array.  The  reflections  which  can  be  seen  in 
Figure  12  as  extra  pips  to  the  sides  of  the  main  lobes,  come  from  aluminized 
mylar  baffles  placed  between  the  filter  elements  above  each  detector 
in  the  array.  The  problem  is  the  greatest  for  the  2.0  -  2.6  pm  center 
element  which  has  these  aluminum  baffles  on  both  sides.  The  result  is 
to  make  each  round  detector  appear  as  two  or  three  detectors  placed  side 
by  side.  Thus,  the  "dip"  that  is  observed  in  the  figures  is  an  optical 
problem  and  its  magnitude  is  a  function  of  the  strength  of  the  side-lobe 
reflection. 

The  data  were  sampled  in  the  cross-track  dimension  at  a  constant 
rate  supporting  digital  reconstruction  of  up  to  .25  cycles/mr.  In  the 
in-track  dimension,  several  sample  frequencies  were  employed  depending 
upon  altitude.  These  frequencies  are  shown  below: 


Altitude 

800 

1600 

3200 

6400 


ln-Trace  Sampling  Frequency  Limit 


.25  cycles/mr 


.50  cycles/mr 
1.0  cycles/mr 
2.0  cycles/mr 


It  is  important  to  know  the  digital  sampling  frequency  with  respect  to 
the  spatial  frequency  response  curves  to  qualitatively  evaluate  the  effect 
of  aliasing.  Aliasing  can  best  be  described  as  a  folding  of  the  spatial 
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frequency  spectrum  at  the  sampling  frequency.  This  phenomenon,  which 
is  a  result  of  the  descrete  time  sampling  of  an  analog  time  continuous 
signal,  causes  spatial  frequencies  higher  than  the  sampling  frequency 
to  appear  to  be  folded  or  "aliased"  to  lower  spatial  frequencies.  If 
these  aliased  signals  have  large  enough  magnitudes,  erroneous  interpretation 
of  the  spatial  information  obtained  from  the  imagery  can  occur. 

The  only  way  to  avoid  aliasing  is  to  make  the  sampling  rate  at 
least  twice  as  high  as  the  highest  expected  spatial  frequency  in  the 
analog  video  data.  Thus,  oversampling  decreases  chance  for  errors,  caused 
by  aliasing  and  the  higher  altitude  data  collected  during  this  program 
can  be  assumed  to  be  aliasing  free  in  the  in-track  dimension.  At  the 
lower  altitudes  in  the  in-track  dimension  and  at  all  altitudes  in  the 
cross-track  dimension,  the  potential  for  some  aliasing  does  exist  as 
can  be  seen  by  examination  of  the  MTF  curves  past  the  .25  cycle/mr 
sampling  frequency  limit.  The  potential  for  aliasing  is  stronger  in 
the  in-track  dimension  (at  least  at  the  lower  altitudes)  than  the  cross 
track.  The  existance  of  an  aliasing  potential  is  a  problem  only  if  the 
scene  contains  significant  structure  at  the  unsampled  frequencies. 
Examination  of  the  reconstructed  imagery  showed  no  noticeable  artifacts 
indicative  of  such  a  condition. 

3.5  RADIOMETRIC  SENSITIVITY 

Measurements  were  made  to  determine  the  radiometric  sensitivity  of 
the  M-7  scanner  in  the  spectral  bands  employed  during  the  NVL  data  collec¬ 
tion  effort.  The  figure  of  merit  commonly  used  to  describe  this  detector/ 
scanner  sensitivity  is  the  noise-equivalent  radiance  difference  (NEAL). 

This  value  is  calculated  by  dividing  the  radiance  difference  producing 
signal  by  the  detector  peak  voltage-to-rms  noise  ratio  of  that  signal 
difference. 

For  the  solar  reflective  bands  in  the  region  from  0.4  to  3.5  pm, 
the  NEAL  is  a  function  of  the  radiance  level,  L.  Hence,  a  NEAL  value  in 
this  spectral  region  really  has  meaning  only  for  a  specified  illumination 
condition.  However,  typical  values  representative  of  the  illumination 
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levels  encountered  during  the  1330  hours  flight  at  A.P.  Hill  were  obtained 
for  the  .8  to  1.1  and  2.0  to  2.6  um  bands.  Shown  in  Table  12  is  the 
digital  reference  source  data  obtained  from  the  fifth  pass  during  the 
1330  hour  mission  for  the  above  two  bands.  For  a  reference  lamp  radiance 
as  shown,  the  NEAL  was  calculated  for  each  band  and  is  shown  in  the 
table. 

For  the  thermal  emissive  region  from  4.0  to  14  um,  the  sensitivity 
of  the  various  bands  can  be  more  usefully  presented  as  noise-equivalent 
temperature  differences  (NEAT).  However,  it  is  important  to  realize 
that  the  rationship  of  radiance  to  temperature  is  non-linear;  hence  a 
given  increment  in  temperature  will  result  in  various  band  radiance  in¬ 
crements  depending  on  the  mean  temperature  and  band  involved.  Shown  in 
Figures  19,  20  and  21  are  the  radiance  versus  temperature  curves  for  the 

2.9  to  3.9,  4.3  to  5.3,  and  8.4  to  13.2  um  bands.  In  Figure  19  it  can  be 

2 

seen  that  a  .02  um/cm  -ster  increment  in  radiance  centered  around  31°6 
results  in  a  AT  equal  to  16°C  while  the  same  radiance  increment  centered 
around  64.5°C  results  in  a  temperature  increment  equal  to  5°C.  Hence,  a 

2 

given  thermal  band  RMS  noise  figure  expressed  in  radiance  units  of  um/cm  /sr, 
will  result  in  different  NEAT's  for  the  system  depending  on  the  mean 
temperature  of  the  measurement.  For  all  of  the  NEAT  calculations  obtained 
in  the  following  tests,  the  mean  temperature  was  assumed  to  be  approxi¬ 
mately  10°C  to  correspond  with  the  surface  scene  temperatures  observed 
during  data  collection. 

The  tests  were  made  using  blackbody  reference  plates  internal  to  the 
M-7's  scanner  housing.  There  are  three  of  these  blackbody  plates  (two 
temperature-controlled,  one  ambient)  which  are  viewed  sequentially  during 
each  scan.  The  purpose  of  these  plates  is  to  provide  known  sources  of 
emitted  radiation  by  means  of  which  the  temperatures  of  objects  viewed 
on  the  ground  by  the  scanner  can  be  calculated.  In  Table  13  are  shown 
the  digital  values  for  the  blackbody  reference  plates  from  the  1330  hour 
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FIGURE  19.  TEMPERATURE  VS  RADIANCE  AX=2.85  -  3.95jjm 


-  13.74m 


flight  and  the  fifth  pass.  Nominal  plate  temperatures  for  this  pass 

* 

were  2.12  C  for  the  cold  plate  and  10.37°C  for  the  ambient  plate  . 

Using  Figures  19,  20,  and  21  and  the  S/N  ratios  shown  in  Table  13, 
the  NEAT  values  for  the  three  detector  bands  at  10°C  are: 

BAND  NEAT  (10°C) 


2.9  - 

3.9 

8. 70°C 

4.3  - 

5.3 

0. 30°C 

8.4  - 

13.7 

0. 27°C 

The  NEAT  of  the  2.9  -  3.9  band  should  be  worse  than  the  other  two  thermal 
bands  since  there  is  much  less  energy  emitted  in  this  band  at  10° C. 
However,  the  difference  observed  in  these  tests  was  larger  than  expected. 
Even  under  daylight  conditions,  where  the  detector  was  receiving  signi¬ 
ficant  solar  reflected  energy  in  addition  to  thermal  emissive  energy 
the  overall  image  quality  was  much  poorer  than  anticipated.  The  con¬ 
clusion  to  be  drawn  is  that  this  element  of  the  array  was  not  performing 
up  to  state-of-the-art  potential,  and  the  resultant  data  are  probably 
of  little  use. 

The  above  discussion  describes  the  nominal  sensitivity  of  the  data 
collected.  Techniques  to  extract  quantitative  thermal  radiometric  infor¬ 
mation  from  any  of  the  data  images  are  given  in  Appendix  C. 


Due  to  digital  hardware  limitations,  the  hot  plate  was  not 
reformatted  onto  CCTs.  Since  only  two  calibrated  thermal  reference 
plates  are  in  fact  required,  this  presented  no  calibration  problems. 
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4.0 

MULTI SPECTRAL  TARGET  CUEING  STUDY 

The  major  activity  of  this  program  was  the  collection,  preparation 
and  documentation  of  a  unique  multitemporal  and  multispectral  data  base 
for  a  tactical  military  scenario.  That  activity  was  successfully  accom¬ 
plished.  However,  the  second  goal  of  providing  significant  initial 
evaluation  of  multispectral  processing  concepts  was  only  partially 
achieved  because  of  unforeseen  diversion  of  resources  to  accomodate  technical 
problems  arising  in  the  data  base  generation  effort.  This  section  describes 
the  analyses  conducted  and  the  results  achieved. 

4.1  GENERAL  CONSIDERATIONS 

The  multispectral  data  base  collected  for  this  program  can  be  considered, 
in  its  fullest  generality,  to  be  an  eight  dimensional  data  set  composed 
of  five  spectral  dimensions,  two  spatial  dimensions,  and  one  temporal 
dimension.  A  full  understanding  of  the  total  information  content  in  this 
data  would  require  a  unified  treatment  of  all  dimensions,  including  second 
order  statistical  interactions. 

Information  in  the  temporal  dimension  is  seldom  considered  in  tactical 
military  applications  of  remote  sensing,  the  obvious  exception  being  in  change 
detection  concepts.  While  temporal  trends  in  spatial/spectral  data  charac¬ 
teristics  may  well  provide  definitive  signature  opportunities,  the  perceived 
need  to  support  real-time  operations  negates  the  potential  utility  of  such 
concepts.  Rather,  time  is  treated  as  a  discrete  ancillary  condition  and 
not  as  a  true  dimension  of  the  information  in  the  data. 

The  nature  of  the  multispectral  scanner  image  provides  two  independent 
spatial  dimensions  to  the  data  set.  The  crudest  use  of  such  data  involves  the 
simple  concept  of  a  pulse  length  in  either  dimension.  More  complex  spatial 
features  such  as  area,  shape  or  texture  can  be  derived.  Use  of  such  derived 
features  tends  to  recognize  the  fact  that,  while  the  data  may  be  independent 
in  an  x,y  coordinate  system,  the  spatial  information  content  of  the  scene 
may  indeed  be  more  uniquely  perceived  by  use  of  a  derived  feature  or 
alternative  coordinate  context. 
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The  simplest  use  of  the  five  spectral  dimensions  of  the  data  is 
just  to  treat  them  as  five  independent  information  channels  in  processing. 
However,  it  is  well  known  that  correlation  does  exist  between  spectral 
bands  in  the  infrared,  so  that  the  true  dimensionality  of  the  information 
is  probably  less  than  five.  Further, the  principal  components  of  infor¬ 
mation  in  the  five  dimensional  spectral  space  do  not  necessarily  line  up 
with  any  of  the  five  axes  (acquired  bands)  in  that  coordinate  system. 
Finally,  much  of  the  information  content  is  really  not  of  interest  (e.g. 
solar  illumination  level  effects)  so  that  by  using  derived  features 
(e.g.  band  ratios)  one  may  be  able  to  differentially  suppress  that 
information  not  of  use  (i.e.  clutter)  in  a  specific  application. 

A  preliminary  evaluation  of  two  different  approaches  to  multispectral 
target  cueing  was  accomplished.  The  first  approach  involved  classical 
maximum  likelihood  classification  using  in-scene  training.  The  second 
approach  involved  unsupervised  spectral/spatial  clustering  followed  by 
discrimination  based  upon  size.  The  results  of  these  efforts  are  pre¬ 
sented  in  the  following  sections. 

4.2  MAXIMUM  LIKELIHOOD  CLASSIFICATION  WITH  TRAINING 

The  most  commonly  used  procedure  for  generalized  multispectral 
recognition  is  that  of  maximum  likelihood  classification  with  in-scene 
training.  In  this  procedure,  the  spectral  values  of  each  pixel  in  the 
scene  are  compared  to  the  statistical  distribution  functions  defining 
the  characteristics  of  the  candidate  classes  from  which  that  pixel  might 
have  been  drawn.  The  pixel  is  then  assigned  to  that  class  from  which 
it  most  likely  was  drawn.  The  candidate  class  distributions  are  empirically 
defined  using  samples  from  the  scene  which  are  somehow  identified,  usually, 
as  in  this  case,  using  ground  truth.  A  general  description  of  this 
methodology  is  given  in  Appendix  E. 

Detailed  procedures  and  results  of  the  maximum  likelihood  classifi¬ 
cation  are  presented  in  Appendix  F.  The  results  are  summarized  in  Table 
14.  It  is  obvious  from  these  results  that  a  significant  lack  of  exten- 
dability  exists  for  the  target  signature  characteristics.  For  the  daytime 
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TABLE  14.  MAXIMUM  LIKELIHOOD  CLASSIFICATION 


Targets 

Present 

Targets 

Cued 

Time 

Training 

Test 

Training 

Test 

0930 

7 

5 

7 

0 

1330 

6 

6 

6 

0 

2330 

3 

9 

3 

6 

data,  the  only  targets  detected  were  those  from  which  the  signature 
training  statistics  were  established.  At  night,  however,  two-thirds 
of  the  test  set  (six  out  of  nine  targets)  were  detected  in  what  amounted 
to  a  thermal  hot  spot  evaluation. 

The  false  alarm  problems  in  the  maximum  likelihood  classification 
were  not  quantitatively  evaluated.  However,  reference  to  the  classifi¬ 
cation  results  in  the  photos  of  Appendix  F  shows  a  qualitative  trend 
quite  similar  to  that  for  target  detection:  false  alarms  are  a  major 
problem  during  the  daytime,  but.  become  significantly  less  so  at  night. 

Had  the  spatial  dimensions  of  the  data  been  utilized  for  additional 
size  discrimination,  the  false  alarms  encountered  at  night  might  well  have 
been  reduced  to  a  potentially  tolerable,  or  at  least  promising,  level. 
However,  a  spatial  discrimination  would  not  have  led  to  such  promising 
results  on  the  daytime  classification  false  alarms. 

4.3  UNSUPERVISED  SPECTRAL/SPATIAL  CLUSTERING  OF  DERIVED  FEATURES 

The  use  of  any  multispectral  recognition  approach  requiring  training 
with  in-scene  ground  truth  is  of  obviously  limited  operational  utility 
in  a  military  reconnaissance  application.  Because  of  this  fact,  unsuper¬ 
vised  approaches  which  attempt  to  partition  the  scene  according  to  the 
natural  spectral/spatial  structure  of  the  data  have  obvious  appeal.  The 
ultimate  utility  of  such  approaches  depends  upon  the  presence  of  recogni¬ 
zable  structure  to  the  data  which  can  be  associated  with  the  presence 
of  a  target  of  interest. 


A  specific  methodology  for  unsupervised  spectral/spatial  clustering 
was  implemented  in  order  to  preliminarily  assess  the  degree  to  which 
useful  and  recognizable  data  structure  might  be  present  in  the  Camp 
A.P.  Hill  data  base.  The  technology  utilized  in  this  investigation  has 
been  successfully  demonstrated  by  ERIM  and  others  in  civilian  applications 
of  multispectral  remote  sensing. 

The  approach  taken  involves  a  preprocessing  of  the  raw  spectral 
bands  in  order  to  create  derived  spectral  features  with  which  to  work. 

This  preprocessing  allows  for  a  more  intuitively  logical,  physically- 
based  understanding  of  the  nature  of  the  data  variability  present  in  the 
multispectral  information.  A  detailed  description  of  the  methodology 
and  results  of  current  state-of-the-art  spectral/spatial  clustering 
applied  to  the  1330  hours  Camp  A.P.  Hill  data  base  is  presented  in 
Appendix  G. 

i 

The  primary  intent  of  the  unsupervised  spectral/spatial  clustering 
investigation  was  to  ascertain  whether  there  existed  any  coherent 
spectral  structure  to  a  tactical  target  which  could  separate  it  from  its 
immediate  background  so  as  to  produce  a  recognizable  spatial  shape.  The 
results  of  this  investigation,  though  admittedly  limited  in  scope, 
suggest  that  such  coherent  spectral  structure  does  not  exist  in  general. 

It  was  found  that  some  targets  clustered  into  recognizable  spatial 
shapes,  while  others  were  completely  inseparable  from  their  backgrounds. 
However,  in  general,  and  for  armor  in  particular,  it  was  found  that  any 
combination  of  the  physically-based  spectral  features  produced  several 
clusters  for  each  target,  so  that  no  single  recognizable  spatial  entity 
could  be  isolated.  That  is,  within  very  broad  tolerance  limits,  targets 
of  interest  do  not  have  coherent  or  uniform  spectral  features  for  the 
spectral  bands  of  the  data  base  investigated. 


5.0 

CONCLUSIONS  AND  RECOMMENDATIONS 


A  mulcispectral  infrared  data  base  of  tactical  military  vehicles 
has  been  created.  Data  were  gathered  at  Camp  A.P.  Hill,  Va.  during 
typical  springtime  conditions  prior  to  leafing  of  deciduous  trees.  The 
background  constituents  were  typical  of  those  found  in  undeveloped  land 
in  temperate  regions.  The  data  base  has  been  documented,  reformatted 
into  computer  compatible  form,  and  delivered  to  the  Sponsor. 

A  limited  analysis  of  the  data  was  accomplished  in  order  to  provide 
initial  indications  of  the  utility  of  multispectral  target  recognition 
concepts.  The  conclusions  reached  in  this  initial  analysis,  and  the 
recommendations  based  upon  them,  are  presented  in  the  following  sections. 

5.1  CONCLUSIONS 

The  existence  of  a  unique,  passive  infrared  spectral  signature  for 
tactical  military  vehicles  has  not  been  verified.  It  appears  that  an 
extremely  large  amount  of  spectral  variability  exists,  both  between 
targets  and,  more  importantly,  even  within  a  single  target.  Attempts 
to  improve  that  situation  by  performing  spectral  transformations  to  isolate 
more  stable  and  physically  interpretable  spectral  features  have  not  been 
successful . 

The  basic  problem  in  the  passive,  spectral  domain  is  that  there 
appears  to  be  no  unique  and  fundamental  cause  and  effect  relationship 
between  what  a  military  vehicle  is,  and  how  it  appears  spectrally  in  the 
infrared  portion  of  the  electromagnetic  spectrum.  A  modest  caveat  on  this 
statement  is  that  some  degree  of  uniqueness,  at  least  with  respect  to 
natural  backgrounds,  does  exist  for  the  thermal  character  of  an  operating 
vehicle  at  night.  Otherwise,  there  may  be  no  reason  to  expect  unique  spectra 
a  priori,  unless  by  human  error  some  infrared  equivalent  to  a  Signal 
Flare  Orange  paint  in  a  snow  background  is  encountered. 

The  problem  in  the  spatial  domain  is  really  a  reflection  of  the 
extreme  amount  of  spectral  variability  present  within  a  given  vehicle. 
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Even  modestly  sophisticated  approaches  to  finding  coherent  spectral 
relationships  generally  fail  to  isolate  a  spatial  entity  which  is  recog¬ 
nizable  by  classical  shape  recognition  concepts.  The  pattern  recognition 
problem,  after  maximal  use  of  the  spectral  information  present,  still 
appears  to  require  the  kind  of  complex  spatial  structure  synthesizing 
capability  typically  provided  by  the  eye/brain  combination  in  human 
vision. 

5 . 2  RECOMMENDATIONS 

The  generally  negative  conclusions  concerning  the  existence  of  some 
useful  spectral  signature  must  be  treated  as  preliminary.  Further  con¬ 
sideration  of  the  problem  may  yield  to  some  creative  Insight  not  presently 
available.  It  is  recommended  that  any  such  investigations  must  recognize, 
and  emphasize  overcoming,  the  basic  problems  faced: 

•  "Tankness"  does  not  uniquely  translate  to  any 
particular  region  or  feature  of  infrared  spectral 
space 

•Different  spectral  bands  are  not  likely  to  resolve 
the  problem,  for  they  are  generically  measuring 
the  same  attributes 

•Different  attributes,  such  as  polarization,  may 
provide  a  key  for  they  may  be  more  closely  dependent 
upon  what  a  vehicle  is  (e.g.,  specular  surfaces). 

Active  techniques  have  been  explored  in  other 
efforts  [4] 

•If  spectral  features  are  found,  they  are  going  to 
be  subtle,  so  that  one  must  concentrate  on  the  real 
targets  (e.g.,  Soviet  armor  in  typical  field  paint 
and  backgrounds),  and  not  upon  questionable  analogs. 

[4]  M.  Bair,  et  al.,  "Active  Passive  Cueing  Techniques"  (U),  EFIM 
Report  No.  112900-69-F,  Air  Force  Avionics  Laboratory,  Contract 
No.  F33615-75-C-1155  (In  Publication,  June  1979)  (SECRET) 


•  The  motivation  for  searching  out  spectral 
features  should  be  driven  by  the  need  to 
isolate  coherent  spatial  entities  amenable 
to  state-of-the-art  spatial  pattern  recognition 
methodologies . 

The  apparent  lack  of  unique  spectral  signatures,  when  coupled  with 
the  extreme  within-target  spectral  variability  found,  suggests  that  the 
means  for  straightforward  definition  of  characteristic  shapes  necessary 
for  classical  pattern  recognition  techniques  may  not  be  attainable.  It 
is  therefore  recommended  that  emphasis  be  placed  upon  developing  pattern 
recognition  technology  which  holds  the  potential  for  emulating  the 
eye/brain  ability  to  synthesize  gross  spatial  structure  through  recognizing 
the  associations  of  the  individual  elements  of  that  structure.  The  three- 
dimensional  target  classifier  (3DTC)  processing  concept  currently  under 
development  at  ERIM  has  investigated  that  potential  [5]. 


[5]  S.  Sternberg,  et  al.,  "3D  Target  Classifier"  (U),  ERIM  Report 

No.  126200-36-F,  Air  Force  Avionics  Laboratory,  Contract  F33615- 
77-C-1010  (In  Publication,  June  1979)  (CONFIDENTIAL) 


APPENDIX  A 

M-7  OPTICAL-MECHANICAL  SCANNER 

The  M-7  optical-mechanical  scanner  is  an  advanced  type  of  airborne 
sensor  capable  of  collecting  and  recording  multiple  channels  of  image 
data  in  either  digital  or  analog  form  on  magnetic  tape  for  subsequent 
processing  and  analysis.  The  processing  can  take  advantage  of  the 
Information  content  of  recorded  spectral  data  in  the  ultraviolet,  visible, 
and  the  thermal  infrared  for  automatic  terrain  classification  and  iden- 
tification.  The  M-7  scanner  has  successf (illy  acquired  high  quality  data 
for  numerous  applications,  including  the  1971  Corn  Flight  Watch  Experiment, 
varied  Landsat  and  Skylab  Investigations,  and  many  others  since  1971; 
its  predecessor  was  first  flown  in  the  mid-1960's.  The  total  system, 
including  boresight  cameras,  can  be  operated  in  ERIM's  C-47,  C-46  or 
C-7  aircraft. 

The  M-7  scanner  intercepts  terrain  radiation  by  means  of  a  5-inch 
diameter  collection  aperture.  This  radiation  is  redirected  onto  an 
optical  telescope  by  a  rotating  flat  mirror.  The  incoming  radiation 
passes  through  a  system  of  beamsplitters,  dispersing  optics,  and  spectral 
filters  to  a  number  of  radiation  detectors  in  the  focal  plate  of  the 
telescope  which  convert  the  focused  beam  of  radiation  to  a  set  of  electrical 
signals  for  recording  on  magnetic  tape.  Figure  A-l  shows  a  schematic 
drawing  of  the  M-7's  optical  system. 

The  rotating  mirror  causes  the  narrow  Instantaneous  f ield-of-view 
of  the  telescope  to  scan  in  a  plane  perpendicular  to  the  longitudinal 
axis  of  the  aircraft  through  an  opening  in  the  bottom  of  the  aircraft. 

Then  before  making  the  next  ground  scan,  it  scans  radiation  reference 
devices  which  are  internal  to  the  scanner.  By  the  time  the  next  scan 
begins,  the  aircraft  has  moved  forward;  thus,  subsequent  line  scans  are 
built  upon  one  another  to  produce  a  continuous  strip  image  of  the  terrain 
beneath  the  aircraft. 
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Covering  a  wavelength  range  from  0.33  to  14.0  micrometers,  the 
M-7  scanner  can  operate  in  up  to  19  different  spectral  bands  through 
the  use  of  various  detector/filter  combinations.  Up  to  12  spectral 
channels  of  image  data  can  be  recorded  at  any  one  time  on  either  a  high- 
density-digital  or  an  analog  14-track  magnetic  tape  recorder.  The 
other  two  tape  recorder  tracks  are  used  for  housekeeping  data  purposes. 

A  typical  grouping  of  detectors  provides  a  common  spatial  resolution  of 
2.5  milliradians .  Other  detectors  have  instantaneous  f ields-of-view 
(FOV)  ranging  from  2.0  to  20.0  mr. 

The  system  has  a  nominal  thermal  resolution  (NEAT)  of  0.1°C 
and  a  nominal  reflectance  resolution  (NEAT)  of  1  percent. 

The  scanner  views  the  terrain  during  90°  of  its  scan,  providing 
an  external  f ield-of-view  (FOV)  covering  the  range  from  45°  left  to 
45°  right  of  nadir.  The  scanner  can  be  tilted  for  oblique  scanning 
ahead  of  the  aircraft  at  angles  of  up  to  52°  from  nadir.  The  scanner 
imagery  is  roll-stabilized. 

The  system  operates  at  either  of  two  constant  scan  speeds  —  60  or 
100  scans  per  second.  Nominal  electronics  bandwidth  is  tape  recorder 
limited  to  a  range  of  dc  to  150  kHz. 


APPENDIX  B 


AIRBORNE  SYSTEM  PERFORMANCE 

This  section  describes  the  various  system  performance  tests,  both 
optical  and  electrical,  performed  on  the  M-7  scanner.  It  was  taken 
verbatum  from  Reference  (6)  and  includes  a  discussion  of  overall  data 
quality  as  well  as  system  performance  figures  and  general  calibration 
procedures. 

B . 1  RADIATION  CALIBRATION 

B.1.1  Visible  and  Near-IR  Radiance  Calibration  (0.4  -  2.6  pm) 

The  task  of  trying  to  provide  absolute  radiance  calibration  in  the 
visible  and  near-IR  regions  for  any  line  scanner  is  difficult,  to  say 
the  least.  Much  has  been  wirtten  about  the  theory  behind  such  radiation 
calibrations,  so  no  attempt  will  be  made  here  to  go  into  any  great  depth 
on  the  subject.  This  section  provides,  as  background,  a  limited  amount 
of  theory  on  radiance  calibration  and  describes  the  procedure  and  results 
of  the  tests. 

Calibration  Tests.  The  scanner  calibration  lamp  signal,  as  described 
in  Section  4.4,  is  used  basically  as  a  transfer  standard  to  obtain  the 
apparent  spectral  radiance  of  terrain  objects.  Its  main  purpose  is  to 
account  for  scanner  system  responsivity  changes  that,  in  general,  are 
almost  impossible  to  monitor  conveniently.  The  scanner  components 
themselves  are  designed  to  respond  linearly  with  increased  radiance  from 
the  ground  terrain  —  i.e.,  an  increase  in  radiance  gives  a  corresponding 
increase  in  detector  signal.  Hence  the  scanner  detector  voltage  for 
each  channel  can  be  represented  by  the  following  general  equation: 


[6].  P.G.  Hasell,  Jr.,  et  al.,  Michigan  Experimental  Multispectral 

Mapping  System,  A  Description  of  the  M7  Airborne  Sensor  and  Its 
Performance,  Report  No.  190900-10-T,  Environmental  Research  Insti¬ 
tute  of  Michigan,  Ann  Arbor,  Michigan,  January  1974. 
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VT+pU)  "  K  LT(X)  +  Lpa)J 


(B-l) 


where 


lt(x) 


PT  Et(X) 


L  (X) 
P 


V  (X)  =  scanner  signal  voltages  recorded 

p  by  aircraft  recorder 

K  =  a  constant  dependent  on  system 
factors  such  as  responsivity  of 
detectors  and  photomultipliers, 
system  transmission  losses,  gains,  etc. 

Lt(X)  =  radiance  of  target  ,  assuming  a 
Lambertian  surface 

L  (X)  =  path  radiance  produced  as  a  result  of 

p  scattering  of  radiation  by  molecular  * 

and  aerosol  particles  in  the  atmosphere 

Pt  =  object  directional  reflectance 

At  present,  the  radiance  qualtity  L  (X)  is  impractical  to  measure 

P 

simultaneously  with  the  airborne  scanner  imagery.  But  to  make  up  for 
this,  considerable  effort  is  going  into  modeling  the  atmosphere  so  that 
the  appropriate  corrections  can  be  made  to  the  data  during  processing. 

In  terms  of  scanner  radiance  calibration,  however,  the  individual 
terms  of  Eq.  (1)  are  less  important  than  the  total  radiance  received  at 
the  scanner's  aperture.  Hence,  the  apparent  target  radiance  (which 
includes  the  path  radiance  contribution)  is  a  more  appropriate  quantity 
to  measure.  Simplifying  Eq.  (1),  then,  we  get: 


vr+p  ‘  KLi+pa) 

LT+p(X)  *  apparent  radiance  of  target, 

Lt(X)  +  L  (A) 

T  P 


(B-2) 


In  each  scanner  channel. 
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Were  it  not  for  the  fact  that  the  constant  K  in  Eq.  (2)  actually 
varies  somewhat  because  of  changes  in  the  detector  and  its  responsivity, 
the  radiance  calibration  of  a  scanner  would  be  easy.  A  quick  look  at 
a  known  radiance  source  such  as  a  reflectance  standard  illuminated  by 
•calibrated  irradiance  standards  would  suffice.  Since  this  "constant" 
does  vary,  however,  a  means  to  remove  K  from  Eq.  (2)  becomes  necessary. 

The  prime  purpose  of  the  lamp  reference  source,  as  stated  earlier, 
is  to  remote  the  factor  K.  The  is  done  by  dividing  the  voltage  obtained 
from  the  target  by  the  lamp  voltage.  Hence,  if  the  target  voltage  is 
given  by  Eq.  (2)  and  the  lamp  voltage  is  given  by 


v*  -  *VX) 


then  by  dividing  Eq.  (2)  by  Eq.  (3)  we  get 


T+p 


V£LT+p(X> 

Lfl(A) 


(B-3) 


(B-4) 


The  radiance  calibration  is  performed  in  a  similar  manner,  except 
that  the  target,  instead  of  being  a  terrain  object,  is  an  object  of  known 
radiance.  If-  the  calibration  target  is  flame-sprayed  aluminum  (a  well 
measured  and  stable  reflectance  standard)  illuminated  by  quartz-iodine 
lamps,  then  the  voltage  generated  by  the  radiance  standard  is 


where  L  (A) 
A  Jo 


vx)  ■ 


|>VJ>  '•] 


(B-5) 


spectral  irradiance  of  quartz-iodine  lamps 

directional  reflectance  of  flame-sprayed  aluminum 

system  constant  which  has  varied  from  K  in 
previous  equation 


Dividing  Eq.  (5)  by  the  lamp  voltage 


Vfcl(X)  -  K’L^A) 


(B-6) 
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we  get 


v  _  WW» 

A  l  L„(A)ir 


(B-7) 


Substituting  Eq.  (7)  back  into  Eq.  (4)  and  rearranging,  we  have 


lt+p<x> 


(B-8) 


Equation  (8)  gives  the  apparent  radiance  of  the  target  in  terms 

of  known  or  measured  quantities.  It  is  obvious  from  Eq.  (8)  that  the 

accuracy  of  radiance  calibration  depends  upon  two  factors:  first,  how 

well  p  and  E  are  known;  and  second,  L  (A)  remaining  spectrally 
AJl  X 

constant.  A  complete  error  analysis  of  M7  scanner  signal  calibration 
is  planned.  Also,  a  more  detailed  explanation  of  radiance  calibration 
and  uses  of  the  sun  sensor  is  given  in  Ref.  [7]. 


Tests  and  Results.  The  calibration  fixture  used  to  calibrate  the 
M7  scanner  is  shown  in  Figure  B-l.  The  lamps  used  as  irradiance  standards 
are  200W  General  Electric  6.6A/T42/1CL  quartz-iodine  lamps  traceable 
to  similar  NBS  standards.  The  reflectance  panel  was  a  sheet  of  flame- 
sprayed  aluminum  painted  with  3M  white  paint.  To  obtain  maximal 
accuracy,  the  3M  panel  was  not  assumed  to  be  a  Lambertian  reflector  but 
instead  its  bidirectional  reflectance  characteristics  were  measured  and 
actually  integrated  over  the  scanner  view  angles  of  interest. 


The  dimensions  shown  in  Figure  B-l  were  obtained  by  analytical 
methods  to  provide  the  maximum  uniformity  of  illumination  on  the  flame- 
sprayed  aluminum  panel,  and  also  by  the  geometrical  restraints  imposed 
by  the  scanner. 


[7]  Hasell,  P.G.,  and  L.M.  Larsen,  Calibration  of  an  Airborne  Multi- 

spectral  Optical  Sensor,  Report  6400-137-T,  Willow  Run  Laboratories, 
Institute  of  Science  and  Technology,  The  University  of  Michigan, 

Ann  Arbor,  Sept.,  1968. 


82 


FIGURE  B-l.  M7  SCANNER  CALIBRATION  FIXTURE 


Figure  B-2  shows  the  result  of  the  integration  of  the  lamp  irradiance 
and  panel  bidirectional  reflectance  to  give  uniform  radiance  from  the 
panel.  In  this  figure,  the  average  radiance  value  viewed  by  the  scanner 
looking  perpendicular  to  the  panel  was  assigned  to  a  value  of  100%, 

and  contours  of  radiance  differing  by  +1%  were  plotted  over  the  rest  of 

* 

the  panel.  The  area  between  the  two  inner  rings  (shaded  area)  represents 
the  area  actually  viewed  by  the  scanner.  The  clear  center  portion  is 
obscured  by  a  mirror  used  for  collecting  thermal  energy.  As  can  be 
seen,  the  radiance  varies  by  -4%  and  +2%  about  the  average  for  a, typical 
spectrometer  channel.  The  total  radiance  L^CA)  in  E9-  (5)  versus 
wavelength  off  the  flame-sprayed  aluminum  panel  is  shown  in  Figure  B-2. 

The  actual  computation  of  the  scanner  calibration  constant  C 
involves  the  last  two  terms  of  Eq.  (8),  i.e., 


C 


(B-9) 


where  V^,  =  signal  voltage  of  scanner  calibration  lamp 

for  each  channel 

VAJ,  =  signal  voltage  of  aluminum  reference  panel 
for  each  channel 

L^(A)  =  total  radiance  of  aluminum  reference  panel 

viewed  by  scanner  for  each  channel. 

However,  the  wide  range  of  the  calibration  lamp  currents  encountered 
during  data  collection  produces  a  series  of  constants  C  corresponding 
to  several  lamp  currents.  Hence  Eq.  (9)  can  actually  be  written: 

ci<x>  '  [Vs)]  (B-l0> 

where  =  radiance  calibration  constant  for  lamp  current 
i  and  wavelength  A. 

To  obtain  the  radiance  of  a  terrain  object  in  channel  n,  the 
appropriate  C^(A)  for  channel  n  at  lamp  current  i  is  found  and  multiplied 
by  VT/Vt.  i.e. 
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SCAN 

DIRECTION 


X  AXIS 


FIGURE  B-2.  RADIANCE  CONTOURS  ON  M7  SCANNER  CALIBRATION  REFERENCE  PANEL. 
View  normal  to  panel  as  seen  by  scanner. 
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V‘>  ■  vx)(^ 

Several  measurements  of  this  constant  C^(X)  have  shown  it  to  change 
whenever  the  scanner  lamp  housing  is  moved,  as  when  the  scanner  is 
removed  from  the  aircraft  for  periodic  maintenance.  For  this  reason, 
several  (^(A)  curves  are  available,  each  of  which  applies  to  a  certain 
time  of  the  year.  To  include  all  of  them  in  this  report  would  serve 
no  practical  purpose,  so  only  a  representative  set  is  shown  in  Fig.  B-3 
parts  a  through  k.  The  error  flags  on  the  curve  are  based  on  a  standard 
deviation  of  +1  as  calculated  by  the  digital  computer,  and  represent 
the  noise  in  the  measurement,  not  the  error  in  radiance  calibration. 

No  radiance  calibration  exists  for  the  UV  channel  nor  for  any  channel 
in  Detector  Position  IB,  because  the  scanner  lamp  does  not  illuminate 
that  portion  of  the  scan  mirror  viewed  by  the  detectors  in  Position  1. 
However,  a  separate  calibration  lamp  for  visible  detectors  in  this 
position  will  be  implemented  in  the  near  future. 

Thermal  Reference  Source  Calibration  (9.3-11.7  yim)  .  As  described 
in  Section  4.4,  the  reference  sources  for  the  thermal  detectors  in  the 
M7  scanner  consists  of  three  blackbody  plates  viewed  sequentially  during 
each  scan.  Two  of  these  plates  are  temperature-controlled  while  the 
third  is  allowed  to  vary  with  the  ambient  temperature.  The  purpose  of 
these  plates  is  to  provide  known  sources  of  emitted  radiation  by  means 
of  which  the  temperatures  of  objects  viewed  on  the  ground  by  the  scanner 
can  be  calculated.  Unfortunately,  a  known  radiation  source  alone  is  not 
enough  to  ensure  accurate  temperature  determinations,  since  several  other 
factors  must  be  taken  into  account. 

(1)  Few  objects  have  an  emissivity  near  1.0.  This  means 
that  only  apparent  temperatures  can  be  measured  by 
an  airborne  scanner. 

(2)  The  atmosphere  through  which  the  radiation  from  the 
ground  is  received  absorbs  and  re-emits  some  of  this 
radiation.  Hence,  unless  atmospheric  corrections  are 
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POWER  RADIATED  TO  SCANNER  ((iW/cm  -sr-jim)  x  100 


8' 


made,  the  apparent  temperature  measurements  only 
apply  at  the  receiving  aperture  of  the  scanner. 

(3)  Perfect  blackbody  reference  sources  are  difficult  to 
construct .  This  means  that  the  reference  sources 
themselves  can  be  viewed  only  a  sources  of  apparent 
temperature. 

The  remainder  of  this  section  outlines  the  steps  taken  to  minimize 
outside  influences  (such  as  the  atmosphere)  on  the  apparent  temperature 
calculations,  and  steps  taken  to  document  for  future  reference  all 
calibration  tests  conducted  on  the  thermal  plates. 

Temperature  Calibration  of  Thermal  Channel.  Many  models  of  atmos¬ 
pheric  transmission  in  the  thermal  region  from  4.0  pm  to  15.0  pm  have 
been  formulated  for  various  meteorological  conditions  in  recent  years. 
These  permit  reasonable  corrections  for  thermal  data  obtained  in  any 
band  in  this  region.  However,  to  minimize  the  number  of  corrections 
needed  to  obtain  temperature  information,  a  9.3  -  11.7  pm  filter  has 
been  used  on  the  thermal  detectors  to  utilize  the  excellent  atmospheric 
transmission  in  that  region.  For  most  purposes,  negligible  atmospheric 
corrections  are  needed  in  this  band  on  dry  days  at  the  altitudes  flown 
by  the  ERIM  C-47  aircraft  (i.e.,  15,000  ft  and  below).  No  attempt 
was  made  to  calibrate  the  reference  plates  at  any  band  other  than  9.3 
to  11.7  pm. 

The  reference  plates  themselves,  as  described  earlier,  are  copper 
plates  which  are  thermoelectr ically  temperature  controlled  and  spray- 
painted  with  a  flat  black  paint.  Although  the  surface  temperatures 
are  actually  controlled  to  within  0.1°C  of  the  desired  temperature  by 
precision  thermistors,  these  plates  do  not  emit  radiation  as  an  ideal 
blackbody;  this  is  primarily  because  of  the  emissivity  of  the  flat  black 
paint.  Measurements  made  on  this  type  of  paint  indicate  an  emissivity 
of  approximately  0.94,  which  means  that  6%  of  the  radiation  coming  off 
a  plate  is  reflected  energy.  Hence,  for  any  plate  temperature  other 
than  the  surrounding  ambient,  the  apparent  plate  temperature  as  seen 
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by  the  scanner  will  be  partly  dependent  on  some  reflected  energy,  and 
will  not,  in  general,  agree  with  the  surface  temperature  as  denoted  by 
the  thermistor. 


Laboratory  tests  were  made  on  the  plates  to  measure  the  difference 
between  surface  temperatures  and  apparent  temperatures.  Strictly 
speaking,  these  measurements  apply  only  to  the  ambient  temperature 
conditions  existing  in  the  laboratory  (approximately  22°-23°C)  at  the 
time  of  the  measurements.  However,  the  apparent  temperature  of  the 
plates  for  any  surface  temperature  can  be  calculated  if  the  ambient 
temperature  in  the  scanner  housing  is  known.  This  is  the  prime  reason 
for  having  an  ambient  plate  as  a  reference. 

Our  measurements  utilized  YSI  (Yellow  Springs  Instrument  Company) 
precision  thermistors  imbedded  in  each  plate  to  control  the  temperature 
through  a  feedback  loop.  These  thermistors  are  calibrated  to  +  0.1°C 
over  the  range  -50°C  to  +100°C.  Calibration  data  giving  the  temperature 
vs.  resistance  relationship  is  shown  in  Table  B-l.  The  radiometric  tem¬ 
peratures  were  obtained  from  a  20°  FOV  radiometer  (Barnes  PRT-5)  calibrated 
against  high  emissivity  (0.99)  conical  blackbodies.  The  result,  shown 
in  Figure  B-4,  is  plotted  as  dial  setting  vs.  plate  temperature,  where 


dial  setting 


thermistor  resistance  (in  ohms) 
100 


The  dashed  curve  is  the  surface  temperature  as  calculated  from  the 
thermistors,  and  the  solid  and  dotted  black  curves  are  the  apparent 
tanperatures  of  the  two  plates.  As  can  be  seen,  the  reference  plates 
are  not  perfect  blackbodies;  some  correction  for  their  emissivity  is 
needed.  The  surface  temperature  curve  does  not  cross  the  apparent  curve 
at  the  ambient  temperature  existing  in  the  lab  (~22°C),  but  instead 
crosses  at  approximately  35°C.  This  is  because  the  radiometer  head, 
being  much  warmer  than  the  surrounding  walls,  increases  the  ambient 
temperature  seen  by  the  plates.  The  emissivity  of  the  plates  as  indicated 
in  Figure  B-4  is  approximately  0.95. 


89 


■HrnwHMmiMHiMwnMNMNlMMMlMMmmmMMaMmi 


! 


TABLE  B-L  TEMPERATURE  VERSUS  RESISTANCE  FOR  YSI  PRECISION  THERMISTOR 


TEMPERATURE  VERSUS  RESISTANCE  -50°  to  +100°C 


TEMP 

RES 

TEMP 

RES 

TEMP 

RES 

TEMP 

RES 

TEMP 

RES 

<°C) 

(Ohms) 

(°C) 

(Ohms) 

(°C) 

(Ohms) 

(°C) 

(Ohms) 

(°C) 

(Ohms) 

-50 

201. IK 

-20 

29.13K 

+10 

5971 

+40 

1598 

+70 

525.4 

49 

187.3K 

19 

27.49K 

11 

5692 

41 

1535 

71 

507.8 

48 

174. 5K 

18 

25.95K 

12 

5427 

42 

1475 

72 

490.9 

47 

162.7K 

17 

24.51K 

13 

5177 

43 

1418 

73 

474.7 

45 

151.7K 

16 

23.16K 

14 

4939 

44 

1363 

74 

459.0 

45 

141.6K 

15 

21.89K 

15 

4714 

45 

1310 

75 

444.0 

44 

132.2K 

14 

20.70K 

16 

4500 

46 

1260 

76 

429.5 

43 

126.5K 

13 

19.58K 

17 

4297 

47 

1212 

77 

415.6 

42 

115.4K 

12 

18.52K 

18 

4105 

48 

1167 

78 

402.2 

41 

107.9K 

11 

17.53K 

19 

3922 

49 

1123 

79 

389.3 

-40 

101. OK 

-10 

16.60K 

+20 

3748 

+50 

1081 

+80 

376.9 

39 

94.48K 

9 

15.72K 

21 

3583 

51 

1040 

81 

364.9 

38 

88.46K 

8 

14.90K 

22 

3426 

52 

1002 

82 

353.4 

37 

82.87K 

7 

14.12K 

23 

3277 

53 

965.0 

83 

342.2 

36 

77.66K 

6 

13.39K 

24 

3135 

54 

929.6 

84 

331.5 

35 

72.81K 

5 

12.70K 

25 

3000 

55 

895.8 

85 

321.2 

34 

68.30K 

4 

12.05K 

26 

2872 

56 

863.3 

86 

311.3 

33 

64.09K 

3 

11.44K 

27 

2750 

57 

832.2 

87 

301.7 

32 

60.17K 

2 

10.86K 

28 

2633 

58 

802.3 

88 

292.4 

31 

56.51K 

-1 

10.31K 

29 

2523 

59 

773.7 

89 

283.5 

-30 

53.10K 

0 

9796 

+30 

2417 

+60 

746.3 

+90 

274.9 

29 

49.91K 

+1 

9310 

31 

2317 

61 

719.9 

91 

266.6 

28 

46.94K 

2 

8851 

32 

2221 

62 

694.7 

92 

258.6 

27 

44.16K 

3 

8417 

33 

2130 

63 

670.4 

93 

250.9 

26 

41.56K 

4 

8006 

34 

2042 

64 

647.1 

94 

243.4 

25 

39.13K 

5 

7618 

35 

1959 

65 

624.7 

95 

235.2 

24 

36.86K 

6 

7252 

36 

1880 

66 

603.3 

96 

229.3 

23 

34.73K 

7 

6905 

37 

1805 

67 

582.6 

97 

222.6 

22 

32.74K 

8 

6576 

38 

1733 

68 

562.8 

98 

216.1 

21 

30.87K 

9 

5265 

39 

1664 

69 

543.7 

99 

200.8 

100 

203.8 

LAMP  RADIANCE  (/xW/cm  -sr./xm) 


FIGURE  B-4.  INTERNAL  LAMP  TRANSFER  STANDARD  CALIBRATION  I  JUNE  TO  30  DE¬ 
CEMBER  1972 


To  a  first  approximation  and,  more  exactly,  for  ambient  housing 
temperatures  around  35°C,  the  curve  shown  in  Figure  B-4  can  be  used  to 
obtain  the  apparent  plate  temperature  calibration  points  for  any  given 
thermistor  control  dial  setting.  Normally,  the  two  temperature-con¬ 
trolled  plates  are  adjusted,  via  oscilloscope  display,  to  bracket  the 
range  of  terrain  signals  viewed  by  the  scanner.  Terrain  temperatures 
are  then  extracted  later  by  interpolation  between  the  apparent  plate 
temperatures.  Care  must  be  taken,  however,  where  there  are  plate 
temperature  differences  of  more  than  20°C  because  although  the  interpola¬ 
tion  used  is  a  linear  procedure,  the  temperature  vs.  radiance  relation¬ 
ship  is  not  linear  in  the  thermal  region.  Figure  B-5  shows  this  relation 
ship  at  X  =  10.5  ym.  As  can  be  seen,  a  linear  approximation  is  possible 
over  a  range  of  20°C,  but  larger  errors  would  appear  at  ranges  greater 
than  20°C. 

For  scanner  housing  temperatures  other  than  35°C,  a  correction  to 
the  apparent  temperatures  shown  in  Figure  B-4  can  be  made.  However, 
for  most  temperatures  in  the  housing,  the  error  in  not  making  this 
correction  is  most  likely  to  be  not  more  than  0.5°C  to  1.0°C.  The 
housing's  actual  ambient  temperature  is  calculated  from  the  ambient 
plate,  which  has  imbedded  in  it  a  YSI  precision  thermistor  identical 
to  the  ones  in  the  controlled  plates.  Because  the  controlled  plates 
are  tilted  at  an  angle  of  45°  to  nadir,  the  ambient  temperature  —  as 
determined  from  the  ambient  plate  —  does  not  really  indicate  the  true 
ambient  radiation  reflected  by  the  controlled  plates.  Some  radiation 
from  the  ground  also  is  reflected  off  the  plates  and  contributes  to  the 
ambient  radiation  coming  from  the  inside  surface  of  the  scanner  housing; 
its  influence  is  small,  however,  and  the  effort  required  to  measure  it 
too  great.  So  far  all  practical  purposes  the  temperatures  determined 
from  the  ambient  plate  are  sufficient  to  make  the  necessary  corrections. 

Accuracy  Test  of  Plate  Calibration  Temperature.  As  a  check  on  the 


plate  calibration  performed  in  the  laboratory,  a  test  utilizing  hot 
and  cold  water  baths  was  devised  to  measure  the  scanner's  overall  tem¬ 
perature  sensitivity  and  accuracy.  Separate  plastic  ice  chests  containing 


FIGURE  B- 5. CONTROL  PANEL  SETTING  VERSUS  APPARENT  TEMPERATURE  (JANUARY  1972) 


heated  and  cooled  water  were  placed  side  by  side  underneath  the  aircraft. 
The  temperatures  of  these  water  baths  were  monitored  constantly  by  both 
a  thermometer  and  a  PRT-5  radiometer.  With  the  scanner  running,  water 
temperatures  were  also  determined  from  the  reference  plates.  These 
three  results  were  compared  for  closeness. 


Figure  B-6  shows  an  oscilloscope  trace  obtained  in  one  of  these 
tests.  The  important  parameters  and  results  of  this  temperature  test 
are  shown  in  Table  B-2. 

As  can  be  seen,  the  agreement  between  calculated  and  measured 
temperatures  is  excellent.  Note  that  no  corrections  were  made  for  the 
ambient  housing  temperature,  since  it  is  fairly  close  to  the  laboratory 
test  conditions  of  35°C  and  the  atmosphere  influence  is  negligible. 
Under  actual  flying  conditions,  however,  both  of  these  effects  would 
have  to  be  considered  to  obtain  anything  approaching  the  accuracy  shown 
above . 


Noise  Equivalent  Temperature  Measurements.  NEAT  information  was 
also  obtained  from  the  data  gathered  in  the  previously  described  experi¬ 
ment.  Generally  speaking,  NEAT  is  defined  as  the  temperature  difference 
required  by  the  thermal  detector  to  obtain  an  RMS  noise  voltage  to  signal 
voltage  difference  ratio  of  unity.  Experimentally,  this  is  obtained 
according  to  the  following  relationship: 


NEAT  = 


where 


Tl 

vl2 

VT]l 

VN,RMS 


temperature  of  plate  2 
temperature  of  plate  1 
signal  voltage  of  plate  2 
signal  voltage  of  plate  1 

RMS  noise  voltage  of  either  plate  2  or  plate  1 
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FIGURE  B-6,  OSCILLOSCOPE  TRACE  OF  9.3  to  11.7  4m 
BAND  IN  M7  SCANNER  DURING  THERMAL  CALIBRATION 


TABLE  B-2.  RESULTS  OF  THERMAL  CALIBRATION  TESTS  IN  9.3-11.7  nm  BAND 

Cold  plate  ref.  temp.  =  24.1°C* 

Hot  plate  ref.  temp.  =  33.7°C* 

Ambient  plate  ref.  temp.  *  28.5°C* 


Thermometer 

PRT-5 

Scanner 

Hot  bath 

28.3°C 

28.3°C 

28.4°C 

Cold  bath 

27.2°C 

26.9°C 

27.2°C 

Ambient  plate 

— 

28.1°C 

‘Obtained  by  thermistor  measurement. 


The  noise  value  is,  of  course,  dependent  on  the  electronic 

bandwidth  of  the  system.  More  specifically,  it  varies  inversely  with  the 
square  root  of  the  bandwidth  for  a  flat  spectral  density  condition.  For 
the  test  case  depicted  in  Figure  B-7, the  bandwidth  was  preamplifier- 
limited  to  approximately  240  kHz.  However,  in  most  cases  the  system  is 
tape-recorder-limited  to  approximately  96  kHz.  Hence,  the  actual  NEAT 
encountered  in  normal  operation  is  found  (from  Fig.  38)  to  be 


NEAT  = 


9.6°C 


0.015  V  /96 
X  0.80  V  X  /24o 


NEAT  =  0.11°C 

It  should  be  mentioned  that  the  noise  voltage,  Vn5rms«  was  ob¬ 
tained  from  the  peak-to-peak  noise  voltage  by  assuming  a  Gaussian  noise 
distribution  and  a  peak  factor  of  8  dB  (which  corresponds  to  peak/RMS 
ratios  of  2.6). 

Spectral  Calibration.  Any  attempt  to  identify  an  object  by  its 
spectral  characteristics  obviously  requires  an  accurate  description  of 
the  spectral  throughput  of  the  measuring  instrument.  In  the  case  of 
the  M7  scanner,  only  the  12-channel  spectrometer  lends  itself  easily 
to  this  kind  of  measurement.  The  detectors  used  in  the  other  two  posi¬ 
tions,  i.e.,  near-IR  and  thermal,  use  bandpass  filters  which  are  an 
integral  part  of  the  detector.  Once  the  filter-detector  characteristics 
are  determined  (they  are  well-documented  upon  purchase),  there  is  little 
need  to  measure  the  spectral  characteristics  unless,  of  course,  some 
type  of  trouble  is  suspected,  such  as  light  leakage,  etc.  The  12-channel 
spectrometer,  on  the  other  hand,  uses  light  pipes  on  the  output  of  a 
dispersing  prism  to  obtain  the  required  bandpasses,  and  is  considerably 
more  subject  to  vibrations,  misalignment,  etc.  Primarily  for  this  reason, 
considerable  effort  is  made  to  maintain  an  up-to-date  spectral  calibration. 

The  spectral  calibration  of  the  12-channel  spectrometer  was  accom¬ 
plished  by  mounting  the  whole  scanner  assembly,  including  spectrometer, 
on  a  collimator  (Figure  B-8).  The  collimator  consists  of  a  large,  20- in. 
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FIGURE  B-7.  COLLIMATOR  FOR  SPECTRAL  CALIBRATION  OF  THE  12-CHANNEL  SPECTROMETER 
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diameter  parabolic  mirror  with  a  100-in.  focal  length.  The  exit  slit 
of  a  monochromator,  serving  as  a  source  of  known  wavelengths,  sits  at 
the  focus  of  the  spherical  mirror.  The  collimated  beam  is  then  directed 
to  illuminate  the  aperture  of  the  scanner.  We  found  that  the  image  of 
the  monochomator  exit  slit  filled  the  entire  scanner  aperture.  But 
because  this  image  did  not  fill  the  instantaneous  f ield-of-view,  calibra¬ 
tion  became  necessary  with  the  scanner  operating  so  that  the  collimated 
beam  was  scanned  across  the  f ield-of-view  which  was  thus  filled  on  the 
average. 

For  this  calibration,  the  entrance  slit  of  the  monochromator 
was  illuminated  by  the  Image  of  a  quartz-iodine  lamp.  The  wavelength 
scale  of  the  monochromator  itself  had  been  calibrated  using  mercury, 
potassium,  sodium,  cadmium,  and  cesium  spectral  lamps,  a  helium-neon 
laser,  and  various  narrowband  filters. 

The  relative  spectral  response  of  the  12-channel  spectrometer  is 
shown  in  Figure  B-8.  Unless  otherwise  noted,  all  wavelength  bands  are 
described  by  their  10%  response  points. 

Spatial  Resolution  and  Registration.  Spatial  resolution  measure¬ 
ments  were  performed  on  the  M7  scanner  using  the  collimator  described 
in  Section  4.2  and  pictured  in  Fig.  17.  Purposes  of  these  measurements 
were  not  only  to  verify  the  spatial  resolution  against  its  calculated 
value  but  also  to  provide  an  overall  check  on  the  alignment  of  the  com¬ 
ponents  in  the  optical  system,  including  registration  in  all  detector 
positions.  The  targets  placed  near  the  focal  plane  of  the  collimator 
were  blackened  panels,  each  with  three  rectangular  holes  near  its  center. 

A  number  of  such  targets  were  constructed  with  the  widths  of  the 
holes  and  spacing  between  them  varied  to  provide  0.5,  1,  2,  4,  10  and 
20  mrad  targets  for  imaging  in  the  detector  plane  of  the  scanner.  For 
a  100- in.  focal  length  parabolic  mirror  and  a  2-mrad  target,  calculations 
show  that  a  hole  width  and  spacing  of  0.2  in.  is  needed.  All  the  other 
targets  were  dimensioned  proportionately. 
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APPENDIX  C 

TEMPERATURE  CALIBRATION  OF  THERMAL  SCANNER  DATA 

A  thermal  scanner  responds  to  radiance  differences  between  objects 
in  the  scene  being  mapped.  The  intervening  atmosphere  and  optical/mechanical 
components  in  the  scanner  itself,  however,  alter  the  signals  received 
in  a  manner  which  makes  absolute  temperature  calibration  difficult  with¬ 
out  the  aid  of  reference  sources.  Procedures  for  obtaining  ground  and 
apparent  target  temperatures  using  reference  sources  will  be  developed 
in  this  appendix. 

A.  Basic  Radiation  Principles 

The  raw  electrical  signals  generated  by  a  thermal  detector  mounted 
in  an  airborne  scanner  are  linearly  proportional  to  the  radiance  (not 
temperature)  emitted  by  the  observed  target  -  modified,  of  course  by 
an  intervening  atmosphere.  The  basic  expression  which  relates  the 
radiance  received  by  the  scanner  to  the  target  radiance  is  given  in 
Equation  (C-l) . 


L  (A)  _  T  (A)  L  (A)  +  L  (A) 

scanner  atmos  target  path 


(C-l) 


where 


scanner 


L 

L 

L 


(A) 

atmos 

(A) 

target 

(A) 

path 


apparent  radiance  at  scanner 
entrance  aperture 

atmospheric  transmission 
radiance  emitted  by  ground  target 
radiance  emitted  by  intervening  atmosphere 


The  radiance  emitted  by  the  target,  L  ^a^get»  depends  on  its 
temperature  and  emissivity  as  shown  in  Equation  (C-2). 


L 


(A) 

target 


T 

(A)  +  (1  -  e)  Lb®  (A) 


(C-2) 
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where 


e 


target  emissivity 


Lgu  ■  radiance  of  a  blackbody  at  target  temperature  TT 
T 

Lgg  O)  =  radiance  of  blackbody  at  ambient  temperature  T 

tT 

Combining  equations  (C-l)  and  (C-2)  and  solving  for  Lgg  (A)  we  get: 


TT 

lbb  <»> 


(X) 

scanner 


L  (X)  -  T 

path  atmos 


{a- 


e) 


T 

L  a 
BB 


<t) 


€  T 


(C-3) 


atmos 


If  all  the  terms  on  the  ri|ht  side  of  Equation  (C-3)  are  known  (thus 
enabling  us  to  determine  LBB  (A)),  then  the  target  temperature  TT  can 
be  found  from  Planck* s  radiation  law. 

T  C 

(A)  "  ~  C  /(AT  )  (C_4) 

X*  (C  *  -  1) 

where 


C  *  2h°2 

ci  "«T 

C2  k 

A  ■  wavelength 

T^,  *  temperature 

h  »  Planck's  constant 

c  •  velocity  of  light 

£1  ■  1  steradlan 

o 

k  -  Boltzman  constant 
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Although  Equation  (C-3)  gives  an  explicit  solution  for  determining 
the  temperature  of  any  target  viewed  by  the  scanner,  in  practice  it  is 
not  easy  to  solve  since  several  of  the  terms  on  the  right  side  of  Equa¬ 
tion  (C-3)  are  extremely  difficult  to  measure  accurately  without  elaborate 

Ta 

ground  teams  and  equipment.  In  most  cases,  estimates  of  (X)  and  c 
are  all  that  is  available  for  ground  targets  while  the  atmospheric  para¬ 
meters  Tatmos  a°d  usually  must  be  determined  from  atmospheric 

models.  The  accuracy  of  the  temperature  calculation  depends,  of  course, 
on  how  many  of  these  parameters  are  available  and  how  accurately  this 
information  is  known.  For  rough  calculations,  an  ideal  atmosphere  and 
target  can  be  assumed  in  which  case  Equation  (C-3)  reduces  to: 


(X) 


scanner 


hi  (A) 


(C-5) 


where  e  and  t  _  have  been  assumed  to  be  1.0  and  L  (X)  *  0. 
atmos  .  .  path 

The  final  unknown,  then,  is  L 

scanner 

B.  Calculation  of  Radiance  at  Entrance  Aperture  of  Scanner 

The  radiance  received  by  the  scanner's  collecting  aperture, 

L  .is  related  to  the  scanner  output  signal  by  the  relationship 

scanner 

given  in  Equation  (C-6) . 


V 


target 


L  (X) 
scanner 


(02)(1X)<»c)(R)(t)  +voff9et 


(C-6) 


where  L  ^  *  radiance  received  by  scanner 

scanner 

ft  *  instantaneous  field  of  view  (IFOV)  of  scanner 
AX  .  wavelength  bandpass 
Ac  -  area  of  collecting  aperture 

R  ■  responsivity  of  detection  (volts/watt) 
t  -  transmission  of  scanner  optics 
Voffset  *  volta8e  offset  in  system  due  to  electronics* 

*Due  to  the  preamplifiers  and  pos tamp lif iers,  power  supplies,  etc.  in 
the  scanner  electronics,  zero  radiance  input  usually  will  not  give  zero 
volts  output. 
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Rewriting  Equation  (C-6)  by  combining  all  the  constants,  we  get: 


V 


target 


k  L 


<*>  +V 

scanner  offset 


(C-7) 


Normally,  the  quantities  8,  AA,  and  can  be  measured  quite  accurately. 
However,  the  detector  resonslvity  R  and  scanner  transmission  r  are 
difficult  to  measure  and  both  can  (and  do)  change  with  time.  For  this 
reason,  blackbody  calibration  plates  are  inserted  into  the  scanner's 
FOV  to  provide  known  radiance  sources  from  which  the  values  of  k  and 
^offset  *n  E9uat*on  (C-7)  can  be  determined. 

In  truth,  perfect  blackbodies  (e  ■  1.0)  are  almost  impossible  to 
construct  and  what  really  are  constructed  are  greybodles  with  emissivi- 
ties  near  1.0.  The  radiance  from  greybody  has  already  been  given  by 
Equation  (C-2) . 


L 


(A) 

target 


elBB  (X)  +  «  -  *>  Si  <X> 


(C-2) 


Hence,  if  the  reference  plat  emlssivities  and  scanner  housing  temperature 
are  known,  then  by  using  Equation  (C-2)  and  two  known  sources  of  radia¬ 
tion,  Equation  (C-7)  can  be  solved  for  k  and  V  ,,  ...  The  solutions  for 

offset 

k  and  VQffset  (obtained  using  two  equations  and  two  unknowns)  is  given 
in  Equations  (C-8)  and  (C-9). 


plate 


-  V 

1  plate 


n 


S 


ti 

ld;  + 

BB 


(1  - 


el) 


BB 


T2 

e2  LBB  + 


(1  - 


e2) 


BB 


} 


(C-8) 


offset  plate  1 


IV  -  V 

^  plate  1  plate  2 


K 


)(ci 


si  ♦ (i  -  ci>  si, 


si +  <i  -  ei>  si 


H 


T2  * 

s  si +  (i  -  e2>  si 


! 


(C-9) 
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If  ■  e 2»  then  Equations  (C-8)  and  (C-9)  reduce  to 


Substituting  Equations  (C-10)  and  (C-ll)  back  Into  Equation  (C-7),  we 
can  solve  for  the  radiance  at  the  scanner  entrance  aperture  in  terms 
of  scanner  voltages  and  blackbody  radiances. 


L  (X) 
scanner 


where  e 

V 

scanner 

V 

plate  1 

V 

plate  2 


plate  emissivity 

output  voltage  of  scanner  viewing  collecting  aperture 
output  voltage  of  scanner  viewing  reference  plate  1 
output  voltage  of  scanner  viewing  reference  plate  2 

radiance  from  blackbody  at  temperature 

radiance  from  blackbody  at  temperature  T ^ 

radiance  from  blackbody  at  scanner  housing 
temperature 
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Equation  (C-12)  can  be  used  in  Equation  (C-l)  or  Equation  (C-5)  to 
obtain  the  desired  ground  target  temperatures.  For  the  M7  scanner, 
the  emissivity  of  the  calibration  plates  is  shown  below: 

£emissivity 

The  following  tables  show  the  reference  plate  temperatures  for 
all  of  the  NVL  Images. 
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MULTISPECTRA L  SCANNER  CALIBRATION  LOG  FOR  MISSION 


Site. 


Camp  A.P .  Hill,  Va. 


TABLE  C-l 


Flight  No.  NVL  #1 

Flight  Date  U2S/7H 
Time  2330  EST 


Scanner  No.  1 

Line 

Run 

Time,  G.m.t., 
hr:min:sec 

Black  Body 
Reference 
Temp..  °C 

— 

Radiation  Reference 
Lamp  Current,  Amps 

Ambient 

Blackbody 

Temp  “C 

Plate  <P3 

Start 

Stop 

Plate  1 

Plate  2 

Lamp  1 

Lamp  2 

1 

ft 

0522 

0523 

4.02 

1  1 

mm 

3.  a 

15,48 

2 

0534 

0535 

4.02 

9.90 

■S 

3.0 

15.87 

3 

0541 

0542 

4.02 

9.90  : 

mm 

3.0 

17.46 

a 

0547 

0548 

4.02 

9.90 

warn 

3.0 

19.25 

5 

1 

0552 

4.02 

9.90 

mm 

3.0 

20.54 

6 

0555 

0556 

4.02 

9.90 

mm 

3.0 

21.32 

■ 

B 

0559 

0600 

4.02 

9.90 

warn 

3.0 

21.70 

■ 

8 

0608 

0609 

4.02 

9.90 

mm 

3.0 

21.19 

■ 

9 

0614 

0615 

4.02 

9.90 

mm 

3.0 

21.10 

■ 

10 

0620 

0621 

4.02 

9.90 

mm 

3.0 

22.06 

■ 

fl 

■ 

B 

■ 

m 

■ 

B 

■ 

B 

■ 

B 

■ 

B 

■ 

B 

MULTISPECTRAL  SCANNER  CALIBRATION  LOG  FOR  MISSION 

Camp  A.P.  Hill,  Va.  NVL 


TABLE  C-2 


Flight  No. 
Flight  Date 


NVL  #2 
3/29/78 

1830  EST 


Scanner  No.  1 


Run 

Time,  G.m.t., 
hr:min:scc 

Black  Body 
Reference 
Temp..  °C 

Start 

Stop 

Plate  l 

Plate  2 

B 

0009 

0010 

6.15 

11.35 

2 

0016 

0017 

6.15 

11.35 

3 

0022 

0023 

6.15 

11.35 

B 

0027 

0028 

6.15 

11.35 

5 

0032 

0033 

6.15 

11.35 

6 

0037 

0038 

6.15 

11.35 

B 

0042 

0043 

6.15 

11.35 

8 

0050 

0051 

6.15 

11.35 

9 

0057 

0058 

6.15 

11.35 

10 

0102 

0103 

6.15 

11.35 

Radiation  Reference 
Lamp  Current,  Amps 


Ambient 

Blackbody 


Plate  3 


11.78 


14.94 


16.29 


17.67 


18.69 


18.33 


17.97 


18.45 


MULTISPECTRAL  SCANNER  CALIBRATION  LOG  FOR  MISSION 


Camp  A. P.  Htll,  Va. 


TABLE  C-3 


Flight  No. 


NVL  "3 


Flight  Date  3/30/78 


0930  EST 


Scanner  No.  1 


Time,  G.m.t., 
hr:min:sec 


Black  Body 
Reference 


Radiation  Reference  , 

Lamp  Current,  Amps  lBlackbody 


Temp..  C 


Plate  l  I  Plate  2  I  Lamp  1 


Lamp  2  I  Plate  3 


MULTISPECTRA L  SCANNER  CALIBRATION  LOG  FOR  MISSION 


Camp  A.P.  Hill,  Va. 
TABLE  C-4 


NVL  04 

Flight  No.  _ 

Flight  Date  _  3/30/78 


Time 


1330  EST 


Scanner  No.  I 


Time,  G.m.t., 
hr:min:sec 


Start 


1838 


1904 


1909 


1913 


1918 


1922 


1926 


Stop 


1859 


1905 


1910 


1914 


1919 


1923 


1927 


1932 


1937 


1941 


Black  Body 
Reference 
Temp..  °C 

Plate  l 

Plate  2 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

2.12 

23.54 

Radiation  Reference 
Lamp  Current,  Amps 


Lamp  2  I  Plate  3 


6.07 


6.07 


6.91 


8.45 


10.37 


11.04 


11.78 


12.05 


11.63 
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APPENDIX  D 

DOCUMENTATION  OF  COMPUTER  COMPATIBLE  MAGNETIC  TAPE  FORMAT 
ROUTINES:  IMSS;  OMSS 

VERSION:  1.0 

DATE:  May  23,  1977 

PROGRAMMER:  J.  Morgenstern 

LANGUAGE:  XTRAN 

SYSTEM:  MTS 

PURPOSE:  Input  and  Output  FSR  to  support  multivariate  data  files 

on  both  mag  tape  and  disc  where  the  number  of  bits  per 
data  value  is  a  multiple  of  8:  8  bits,  16  bits  and 

32  bits/value  are  available.  It  also  supports  variable 
length  ancillary  data  fields  at  the  beginning  of  each 
scan  line,  as  well  as  packing  the  data  by  pixels  or 
by  channels. 

INTERFACE:  11LINE/QLINE  System 

COMMENTS:  1)  IMSS  can  only  be  used  to  read  data  files  written 

using  OMSS. 

2)  The  ancillary  data  is  stored  in  a  new  common  block: 

COMMON /QANC I L/QHDNUM , QHEAD ( 2  00 ) 

(QHDNUM  is  the  number  of  bytes  of  ancillary  information) 

3)  OMSS,  at  the  present  time,  takes  on  the  responsibility 
for  ascertaining  the  parameters  which  define  the 
manner  in  which  the  data  are  written  out.  The 
routine  uses  a  NAMELIST  1/0  to  request  the  following 
variables: 

Q0BYTE  :  Number  of  bytes  to  use  to  output 

each  data  value  (may  be  1,  2  or  4) 

Q6IIEAD  :  Number  of  bytes  of  ancillary  information 
at  beginning  of  each  scan  line. 

INTER  :  .TRUE.  Write  out  as  channel 
interleaved  data 

.FALSE.  Write  out  ctiannel  by  channel 

TAPE  :  .TRUE.  Output  will  go  on  tape 

•FALSE,  (default)  output  will  go 
to  disc  file. 

If  the  input  data  is  read  using  I8BIT,  then  the  default  values  for 
the  first  three  parameters  will  he  those  used  on  the  input  data.  If  the 
input  data  is  in  some  other  format,  the  defaults  are  1,  0,  and  TRUE, 
respectively. 
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SYSTEM  CONSIDERATIONS 


A.  New  common  blocks:  QANCIL  and  QI6PAR  were  defined.  The  values 
in  QI9PAR  should  eventually  he  incorporated  in  QINCOM. 


OMSS  keeps  track  of  the  maximum  line  length  available  on  the 
output  device  (32767  bytes  to  MTS  disk,  3000  bytes  to  "tape) . 
When  the  output  line  is  longer  than  this  maximum,  it  is  broken 
up  into  several  pieces  and  multiple  records  per  scan  line  are 
written. 


C.  File  Definition  (Format) 


A.  First  record  is  TITLE/HEADER  record  (it  largely  follows  the 
7094  TITLE  BLOCK). 


Bytes 

Name 

1-4 

NSIZE 

5-8 

QREEL 

9-12 

QFILE 

13-16 

DLINE 

17-20 

QNSA 

21-24 

QNSB 

25-28 

QKS 

29-32 

QNSS 

33-36 

QNA 

37-40 

QKP 

41-44 

QNCILAN 

45-48 

QBYTES 

49-52 

QHEAD 

53-56 

QRECS 

Function 


Number  of  bytes  needed  for  input  buffer 
(scan  line  length) 


(for  multi-tape  files) 

File  number  (for  multi-file  tapes) 


Digital  data  set  line  number  of  first 
scan  line  in  file 


Sensor  (permanent)  number  of  first  line 
in  file 


Sensor  (permanent)  number  of  presumed 
last  line  in  file 


Line  skip  factor  for  stored  data 
Number  of  pixels  per  scan  line 
Number  of  first  scan  point 
Point  numbering  increment  of  data 
Number  of  channels  of  data 


Number  of  bytes  per  individual  data 
value 


Number  of  bytes  of  ancillary  information 
at  beginning  of  scan  line  [multivariate 
data  begins  at  byte  (16+QHEAD)+ll 


Number  of  records  used  to  store  one 
scan  line  of  data 


57-60 

QINTER 

Interleaved  flag: 

0  channel  interleaved 

1  data  stored  by  pixels 

61-64 

QUN1TS 

Number  of  individual  data  values 
stored  in  each  record 

65-68 

DGEN 

Data  set  generation  number 

69-76 

DATE 

Date  of  data  set 

77-140 

SPARE 

(Extra  locations  for  future  use) 

141-260 

QTITLE1 

120  character  title 

481-600 

QL1ST 

120  character  list  of  modules  which 
have  modified  this  data  (6  characters 
per  module  name) 

601-720 

TI1LE2 

Auxiliary  title  block 

721-1440 

TANCIL 

Free-form  ancillary  data  block  to  be 
further  defined  later 

Scan  Line 

Format 

The  remaining  records  in  the  file  are  scan  lines,  with  each 
scan  line  written  in  QREC  consecutive  records.  The  format 
of  any  record  is: 

Bytes 

1-4  Digital  (data  set)  line  number 

5-8  Sensor  (permanent)  line  number 

9-12  Record  number  for  this 

13-16  Blank 

[if  this  is  record  1,  then  bytes  17  through  (QHEAD+16)  ancillary  data3 
The  data  follows  from  here. 

H  / - 

•  •  • 

n\ - 


QINTER  =  0  means: 


Chan 

1 

Chan 

2 

Chan 

3 

Chan 

1 

Chan 

2 

Pixel 

1 

Pixel  2 

QINTER  =  1  is  the  case: 


Pixel 

Pixel 

Pixel 

— 7  ' 

Pixel 

Pixel 

Pixel 

/ 

1 

2 

3 

•  • 

-J  L- 

•  •  • 

QNSS 

1 

2 

/ 

•  •  • 

• 

Channel  1 


Channel  .2  Channel  3 
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When  the  data  must  be  written  out  using  several  physical  records 
per  scan  line  (a  logical  record),  the  following  rules  are  observed; 

a)  All  records  written  contain  the  same  number  of  bytes, 
NSIZE  (which  is  a  multiple  of  180  so  that  the  tapes 
may  be  read  on  any  machine) .  The  records  are  padded 
with  trailing  zeros  after  the  last  data  value  on  each 
record.  The  number  of  data  units  (one  pixel,  one  channel) 
per  record  is  given  in  the  variable  QUNITS  in  the  title 
block. 

b)  The  data  are  broken  between  successive  records  at  a 
boundary  between  data  units:  pixels  (if  the  data  are 
scored  channel  interleaved)  or  channels  (if  the  data 
are  stored  by  channels). 
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APPENDIX  E 


AN  INTRODUCTION  TO  MULTISPECTRAL  CLASSIFICATION 
AND  RECOGNITION  PROCESSING 

Multispectral  classification  essentially  consists  of  deciding  the 
class  of  the  object  being  observed  on  the  basis  of  its  color.  While 
the  human  eye  utilizes  three  dimensions  of  color  —  value,  hue,  and 
intensity  —  for  classification,  other  multispectral  classification 
systems  may  use  as  many  dimensions  as  the  number  of  spectral  bands 
observed.  In  the  following  paragraphs  a  brief  discussion  of  the  theory 
of  multispectral  classification  and  recognition  processing,  as  it  re¬ 
lates  to  the  outputs  of  a  multispectral  scanner,  is  presented.  In 
order  to  simplify  the  explanation  of  the  recognition  processing,  much 
of  the  discussion  will  be  presented  first  in  terms  of  two  spectral 
bands  and  the  extension  to  many  bands  will  be  explained  at  appropriate 
points  in  the  development.  Also,  for  simplification,  we  assume  that 
only  three  materials  are  being  scanned. 

First,  then,  let  us  consider  the  nature  of  two-band  scanner  data 
which  would  be  obtained  from  each  of  three  separate  materials.  Suppose 
that  the  spectral  radiation  intensity  received  at  the  scanner  from  each 
of  the  three  materials.  A,  B,  and  C,  is  as  shown  by  Figures  E-la,  lb, 
and  lc,  respectively.  The  specific  wavelengths  and  A^  indicated  on 
each  of  the  three  figure  parts  correspond  to  the  centers  of  the  respec¬ 
tive  wavelength  bands  covered  by  the  two  bands  of  the  scanner,  and 
therefore,  for  each  material,  the  response  of  a  given  scanner  band  will 
be  proportional  to  the  height  of  the  material's  spectral  curve  at  the 
wavelength  corresponding  to  that  channel.  Thus,  if  u^  is  the  signal 
from  band  1  and  u^  is  the  signal  from  band  2,  the  relative  magnitudes 
of  u^  and  for  each  of  the  three  materials  will  be  as  indicated  in 
Figures  E-la,  lb,  and  lc. 

The  scanner  responses  for  the  three  materials  may  be  presented  in 
a  more  compact  form  by  considering  u^  and  u^  as  the  two  components  of  a 
two-dimensional  vector  and  plotting  the  coordinates  for  each  material 
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FIGURE  E-2.  RESPONSES  OF  TWO  CHANNELS  FOR  MATERIALS  A,  B,  AND  C 


FIGURE  E-3.  RESPONSES  TO  MANY  SAMPLES  OF  MATERIALS  A,  B,  AND  C 
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as  shown  in  Figure  E-2.  The  u^,  plane  shown  will  be  referred  to  as 
signal  space  or  "u"  space.  If  the  scanner  had  three  bands  instead  of 
two,  this  space  would  be  three-dimensional,  with  the  response  of  the 
third  band  corresponding  to  the  third  dimension.  If  the  scanner  had  n 
bands,  the  corresponding  "u"  space  would  be  n-dimensional.  Although  an 
n-dimensional  space  for  n  greater  than  3  is  difficult  to  visualize,  it 
may  be  easily  described  and  handled  mathematically. 

Actually,  because  of  various  statistical  fluctuations  in  the  pro¬ 
perties  of  the  materials  being  scanned,  the  intensity  and  spectral 
properties  of  the  illumination,  the  scanner’s  look  angle  and  other  fac¬ 
tors,  the  plot  of  u^  vs.  u^  will  not  generally  be  distinct  points  as 
indicated  in  Figure  E-2.  Instead,  if  points  for  a  large  number  of 
samples  of  these  materials  are  plotted  in  u  space,  the  points  will  tend 
to  form  three  clusters  as  shown  in  Figure  E-3  with  each  cluster  corres¬ 
ponding  to  one  of  the  three  materials.  A,  B,  or  C.  In  general,  the 
density  of  points  will  be  greater  near  the  center  of  each  cluster  and 
will  become  very  low  near  the  edge.  Also,  the  clusters  will  tend  to  be 
elliptical  rather  than  circular  because  of  correlation  between  changes 
in  u^  and  changes  in  Uj  for  a  given  material.  This  means,  simply,  that 
if  u^  increases  because  of  some  natural  occurrence,  such  as  an  increase 
of  illumination  on  the  scene  being  scanned,  Uj  will  probably  increase 
also. 

The  problem  to  be  solved  by  the  processor  may  be  stated  as  follows: 
"Given  any  sample  point  on  the  u^,  plane,  from  what  type  of  material, 
A,  B,  or  C,  was  the  sample  most  likely  obtained?"  If  the  sample  point 
falls  near  the  centroid  of  one  of  the  clusters  of  points  for  A,  B,  or 
C,  the  decision  is  obvious;  the  material  probably  belongs  to  the  class 
indicated  by  the  group  near  whose  centroid  the  sample  point  is  located. 
Suppose,  however,  that  the  sample  point  is  located  at  e  in  Figure  E-3 
and  thus  does  not  clearly  belong  to  either  A,  B,  or  C.  A  decision  can 
still  be  made  by  comparing  the  densities  (also  called  "likelihoods") 
of  points  from  material  4»  from  material  B,  and  from  material  C  in  the 
neighborhood  of  point  e.  * 
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Assume  that  a  large  strip  of  terrain  has  been  scanned,  and  that 

the  resulting  large  nunfcer  of  sample  points  has  been  plotted  in  the 

ui»  u2  Plane  as  ln  Figure  E-3.  Let  D^Uj.,  u2>  be  the  density  of  sample 

points  from  material  A  as  a  function  of  u  and  u,,  let  D_(u, ,  u.)  be 

1  2  B  1  2 

the  density  of  sample  points  from  material  B,  and  let  D (u. ,  u,)  be  the 
density  of  sample  points  from  material  C.  The  "maximum  likelihood" 
decision  rule  is  to  decide  that  point  e  belongs  to  the  material  with 
the  largest  density  at  e.  The  maximum  likelihood  principle  is  highly 
respected  in  its  own  right.  It  can  also  be  justified  as  a  Bayesian  de¬ 
cision  with  equal  prior  probabilities  and  equal  costs  of  misclassiflcatlon. 

In  a  target-background  type  of  decision  rule,  the  density  of  one 
material,  say  A,  is  compared  with  a  composite  density  of  the  other 
materials,  and  the  rule  is 


decide  on  A  if 


kada 


+  KcDc 


>  1 


where  KA>  Kg,  and  Kc  are  the  relative  frequencies  (prior  probabilities) 
of  A,  B,  and  C,  respectively. 


The  recognition  processor  can  be  implemented  such  that  it  decides 
that  a  sample  belongs  to  a  given  material  (A,  for  example)  if  the  "A" 
likelihood  ratio  for  the  sample  is  greater  than  some  constant,  K. 


Thus,  if 

D.(e) 

- " -  >  K 

Ve)  +  Dc(e)  =  * 


(E-l) 


then  the  processor  decides  that  the  sample  is  material  A.  The  value  of 
the  constant  K  in  equation  E-l  depends  upon  the  estimated  costs  for 
false  alarms  and  missed  detections  and  upon  the  particular  optimal 
decision  criterion  which  is  to  be  employed.  For  recognition  of  various 
crops,  a  value  of  unity  has  been  used  for  K.  In  this  case,  the  identi¬ 
fied  material  must  have  a  point  density  at  point  e  of  at  least  50%  of 
the  total  point  density  for  all  materials  at  e. 
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In  order  to  implement  equation  E-l,  it  is  necessary  to  generate 

mathematical  expressions  which  accurately  describe  D  (u  ,  u.),D_(u  ,  u  ), 

A  1  2.  d  1  2 

and  Dq(u2»  u2^*  is  assumed  that  these  densities  of  points  on  the 
u^,  u2  plane  may  be  represented  in  terms  of  bivariate  Gaussian  probabil¬ 
ity  density  functions  in  this  plane,  if  the  latter  are  fitted  to  the 

measured  statistical  properties  of  actual  data.  Thus,  if  f.(u, ,  u„)  is 

A  1  L 

the  bivariate  Gaussian  probability  density  function  which  represents  the 
relative  density  distribution  of  points  from  material  A  in  the  u^,  u2 
plane  (Fig.  E-3),  the  actual  density  as  a  function  of  u^  and  u2  will  be 
given  by 

W  u2>  “  PCA)^^,  u2)  (E-2) 

where  P(A)  is  the  probability  of  occurrence  of  material  A  in  the  scanned 
area.  Similarly, 

W  u2*  “  u2>  (E-3) 

and 

DC(U1,  u2>  -  P(C)£C<U1,  u2>r  (E-4) 

/ 

! 

Substituting  E-3,  E-4,  and  E-5  into  E-2  gives 
DA(e)  P(A)fA(e) 

DB(e)  +  Dc(e)  =  P(B)fB(e)  +  P(C)fc(e)  ~  K  (E"5) 

if  point  e  is  to  be  classified  as  belong  to  material  A. 

Inequality  E-5  is  the  relationship  which  is  implemented  on  the 
recognition  processor.  Instead  of  working  on  two  bands,  the  processor 
can  accept  information  from  many  channels  and  hence  the  recognition 
space  is  a  multi-dimensional  space  Instead  of  the  two-dimensional  plane 
shown  in  Figure  E-3. 

For  processing  signals  from  an  unknown  area,  the  probabilities  of 
occurrence  of  the  various  materials  (i.e.,  P(A),  P(B),  and  P(C)  in 
eq.  E-5)  must  be  estimated,  since  they  are  not  known.  Experience  has 
shown,  however,  that  these  estimates  can  be  very  crude  without 
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appreciably  affecting  the  processor  operation.  In  fact,  the  processor 
usually  operates  quite  satisfactorily  with  P(A)  =  P(B)  *  P(C).  Equal 
prior  probabilities  were  assumed  for  the  A.  P.  Hill  processing. 

Having  determined  which  of  many  materials  the  unknown  sample  is 
most  like,  the  question  still  remains  whether  the  sample  is  to  be 
classified  as  that  material.  It  is  possible  that  although  the  inequality 
E-5  for  that  material  was  satisfied,  the  magnitude  of  the  probability 
density  function  at  point  "e"  was  extremely  small.  This  indicates  that 
the  sample  at  "e"  might  better  be  classified  as  not  belonging  to  any 
material  whose  probability  distribution  functions  was  included  in  the 
classification.  Therefore,  the  final  classification  or  lack  of  classi¬ 
fication  of  a  sample  is  based  on  its  existence  within  a  fixed  region  of 
the  probability  distribution  function.  This  region  is  very  often  denied 
as  that  region  which  includes  99.9%  of  the  points  in  the  distribution 
function. 
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APPENDIX  F 

SUPERVISED  MULT I S PECTRAL  CLASSIFICATION 

The  supervised  multispectral  classification  was  performed  using 
standard  algorithms  and  equipment  available  in  ERIM's  Earth  Resources 
Data  Center.  A  maximum  likelihood  recognition  algorithm  (see  Appendix  E) 
was  utilized  on  a  DEC  PDP  11/70  computer.  Selected  runs  of  A.  P.  Hill 
data  representing  the  best  spatial  resolution  (1.6  ft)  and  nadir  viewing 
conditions  were  processed.  The  methodology  and  results  of  the  signature 
training  and  subsequent  classification  are  described  in  the  following 
sections. 

F.l  SIGNATURE  TRAINING 

Training  sets  were  selected  by  identifying  the  four  comers  of  a 
quadrilateral  which  defined  each  target  and  background  area  to  be  used 
for 'training.  The  training  areas  used  in  this  study  are  documented  in 
Tables  F-l,  F-2,  and  F-3  for  the  0930,  1330,  and  2330  flight  times, 
respectively.  These  tables  show  the  coordinates  of  the  comers  that 
define  the  training  set  region,  the  name  provided  each  training  set 
and  the  number  of  picture  elements  (pixels)  in  each  training  set.  The 
"X”  coordinate  is  the  pixel  number  location,  whereas  the  "Y"  coordinate 
is  the  scan  line  number  location  within  the  various  data  bases.  The 
vehicle  notation  refers  to  the  target  location  numbers  defined  for  the 
A.  P.  Hill  test  site  rather  than  to  a  specific  target  or  vehicle.  It 
was  originally  planned  that  the  data  gathered  at  1830  be  processed  as 
well.  However,  because  of  limited  target  contrast  for  this  data  base, 
the  targets  could  not  be  dependably  located  and  this  data  set  could  not 
be  included  in  the  processing. 

The  mean  signal  values  and  covariance  statistics  were  established 
for  the  targets  and  backgrounds  selected  for  training.  Tables  F-4, 

F-5,  and  F-6  show  the  statistical  characteristics  of  each  training 
signature  for  the  0930,  1330,  and  2330  hours  data  bases,  respectively. 

For  the  2330  data  only  Channels  4  and  5  (the  two  thermal  channels)  have 
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useful  data  since  there  was  Inadequate  signal  available  In  the  other 
spectral  regions.  For  the  0930  and  1330  data  bases,  all  four  Infrared 
channels  are  available  for  consideration.  The  values  In  these  tables 
are  the  raw  scanner  digital  output  values.  Multispectral  recognition 
using  in-scene  training  does  not  require  a  calibrated  data  base. 

However,  calibration  information  for  these  data  bases  has  been  developed 
and  is  available.* 

Upon  examining  the  signature  statistics  in  Tables  F-4,  F-5,  and 
F-6,  two  interesting  points  are  noted.  First,  the  largest  separations 
begween  the  target  classes  and  background  classes  occur  for  Channels  4 
and  5,  the  thermal  channels.  Also,  the  variation  in  the  signals  attri¬ 
buted  to  the  target  classes  are  much  larger  than  those  associated  with 
the  backgrounds . 

There  are  two  factors  that  may  account  for  this  latter  fact.  The 
first  has  to  do  with  the  difficulty  of  accurately  defining  the  training 
set  area  associated  with  each  target.  It  is  probable  that  considerable 
background  signal  contamination  occurs  within  the  target  training  sets. 
The  second  factor  relates  to  the  fact  that  the  targets  themselves  exhibit 
considerable  variation  in  radiance  over  their  spatial  extent.  This  leads 
to  a  large  amount  of  variance  in  vehicle  training  signals,  at  least  for 
the  fairly  high  resolutions  processed. 

After  the  training  process  had  been  completed,  an  upper  bound  for 
potential  classification  accuracy  was  determined  by  using  the  training 
signatures  to  classify  the  training  data  itself.  Tables  F-7,  F-8,  and 
F-9  present  the  classification  results  for  the  0930,  1330,  and  2330 
data  bases,  respectively.  The  tables  show  the  percent  of  a  specific 
training  set's  pixels  that  are  classified  as  belonging  to  each  signature 
class.  For  example,  in  Table  F-7,  training  set  3  (Trees  2)  is  correctly 


*Beard,  J.  L.  and  R.  Horvath,  "Sensor/Seeker  Subsystems  Modeling," 
Report  No.  112900-64-F,  Environmental  Research  Institute  of  Michigan, 
Ann  Arbor  [in  publication]  (CONFIDENTIAL). 
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classified  98%  of  the  time  as  training  set  3,  1%  of  the  time  it  is 
classified  as  either  training  set  4  (Vehicle  2)  or  training  set  8 
(Vehicle  12).  In  other  words,  the  tables  provide  the  probability  of 
correct  classification  and  false  alarm  rates  within  the  training  sets 
themselves.  These  data  show  that  the  training  signatures  are  fairly 
distinct.  For  this  set  of  data,  the  maximum  confusion  rate  between 
targets  and  background  is  5%  and  occurs  for  the  2330  hours  data.  An 
error  rate  of  5%  represents  one  target  pixel  being  misclassif led  as 
a  background  pixel. 

F.  2  SCENE  CLASSIFICATION 

The  next  step  in  the  processing  was  to  apply  the  training  signa¬ 
tures  from  a  particular  run  to  the  entire  data  base  for  that  same  run 
in  order  to  determine  the  performance  of  this  multispectral  technique 
in  terms  of  both  false  alarms  and  detection.  The  results  of  these 
multispectral  classifications  are  shown  in  Figures  F-l,  F-2,  and  F-3 
for  the  0930,  1330,  and  2330  data  bases,  respectively.  These  figures 
are  photographs  of  the  display  of  the  classification  results. 

In  order  to  interpret  these  photos,  one  needs  the  information 
related  to  their  generation.  The  associated  color  codes  used  for  each 
classification  are  provided  as  a  part  of  Tables  F-l,  F-2,  and  F-3  under 
the  column  labeled  "COLOR".  The  three  digits  shown  under  this  column 
are,  respectively,  the  amount  of  red  (first  column),  the  amount  of 
green  (second  column),  and  the  amount  of  blue  (third  volumn)  used  in 
the  false  color  presentation.  For  example  in  Table  F-l,  the  trailing 
set  labeled  "Vehicle  4"  has  a  color  code  700.  The  display  has  seven 
levels  of  intensity  for  each  color,  thus  700  is  bright  red  in  the  false 
color  display. 

For  the  most  part,  the  vehicle  training  sets  have  been  assigned 
a  bright  red  color,  trees  are  a  shade  of  green,  and  soils  are  white. 

In  the  classification  photos,  pixels  appearing  as  black  are  those  which 
have  been  classified  as  not  belonging  to  any  of  the  training  sets  (i.e., 
they  lay  outside  of  the  99.9%  probability  contours  of  all  signatures). 
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Upon  examining  Figure  F-l,  it  is  seen  that  seven  vehicles  are 
cued.  However,  there  are  a  significant  number  of  false  alarms  within 
the  frame,  especially  within  the  wooded  areas.  This  is  explained  by 
referring  back  to  Table  F-7,  where  it  is  shown  that  a  misclassif ication 
occurred  between  the  categories  Forest,  Vehicle  2,  and  Trees  2.  For 
the  most  part,  the  false  alarms  that  occurred  did  not  have  as  large  a 
spatial  extent  as  the  vehicles  that  were  cued.  Thus,  a  spatial  filter 
could  be  used  to  eliminate  some  of  the  false  alarms.  There  are,  however, 
a  number  of  false  alarms  that  have  spatial  extents  equivalent  to  the 
vehicles  and  as  such  would  not  be  eliminated  using  a  simple  spatial 
filter.  Even  given  such  a  capability  for  false  alarm  suppression, 
however,  a  significant  detection  problem  still  exists  since  only  the 
seven  vehicles  used  for  training  were  identified. 

For  the  1330  hours  data  base  (Fig.  F-2),  six  vehicles  were  used 
for  training.  There  were  no  misclassif ications  between  vehicle  and 
non-vehicle  training  sets  (see  Table  F-8) .  However,  again  in  this  case 
only  the  trained  targets  were  cued  and  in  addition  there  were  more  large 
area  false  alarms  produced  than  are  apparent  in  the  0930  hours  data  base. 

For  the  2330  data  base  (Fig.  F-3)  where  only  3  vehicles  were  used 
for  training,  9  vehicles  were  cued.  False  alarms  tended  to  be  spatially 
small  for  this  data  base.  The  red  lines  that  cross  the  images  are  a 
result  of  line  drop  outs  in  the  original  data  base,  thus,  should  be 
ignored  when  interpreting  the  results.  It  is  also  interesting  to  note 
that  while  the  ambient  M-60  tank  (Position  12)  was  not  used  for  training, 
it  is  one  of  the  targets  detected.  Further,  it  should  be  noted  that 
for  this  data  base  only  the  4.5  to  5.5  and  8  to  14.0  ym  bands  are  being 
used.  Thus,  for  this  data  base,  where  the  targets  tend  to  be  much 
hotter  than  the  background,  a  two-channel  hot  spot  detection  performs 
quite  well. 

The  results  of  this  preliminary  multispectral  cueing  study  can  be 
summarized  as  follows: 
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(1)  The  false  alarm  rate  is  high  for  the  three  data 
bases  examined. 

(2)  For  the  0930  and  2330  hours  data  bases,  a  spatial 
filter  could  significantly  lower  the  false  alarm 
rate. 

(3)  For  the  most  part,  only  those  targets  which  were 
included  within  the  training  set  were  detected 
(exception  was  the  2330  data  base  for  which  9  targets 
were  detected  while  only  3  were  used  for  training). 
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APPENDIX  G 

UNSUPERVISED  SPECTRAL/SPATIAL  CLUSTERING 


Unsupervised  spectral/spatial  clustering  was  performed  using 
standard  methodology  and  algorithms  which  have  been  developed  by  ERIM 
for  earth  resource  and  agricultural  remote  sensing  applications.  The 
basic  approach  consisted  of  three  steps:  1)  transformation  of  the  original 
spectral  channels  into  derived  channels  which  (hopefully)  more  uniquely 
represent  interpretable  and  r  able  information  features  of  the  scene; 

2)  unsupervised  clustering  of  these  features  to  produce  "blobs"  (i.e. 
contiguous  localized  groups  of  pixels  exhibiting  similar  spectral 
characteristics);  and  3)  discrimination  among  blobs  based  upon  their 
target-like  size. 

The  methodology  for  implementing  this  approach  was  intentionally 
constrained  to  use  only  such  information  as  might  be  anticipated  a  priori 
based  upon  physical  principles,  or  which  might  be  derived  from  a  statis¬ 
tical  treatment  of  the  image  data  without  knowledge  of  ground  truth. 
Detailed  methodology  and  results  of  application  to  the  1330  hours  A.P. 

Hill  data  (Run  5)  are  presented  in  the  following  sections. 

G.l  SPECTRAL  FEATURE  EXTRACTION 

Four  spectral  bands  were  utilized  in  this  procedure:  0.86  to  1.04  pm, 
1.9  to  2.55  pm,  4.3  to  5.3  pm,  and  8.4  to  13.7  pm.  The  2.9  to  3.9  pm 
band  was  excluded  because  of  noise  problems  previously  cited  (Section  3.5). 

The  two  near- inf rared  spectral  bands  respond  to  a  combination  of 
solar  illumination  and  object  reflectance  factors.  The  former  factor 
is  highly  correlated  between  these  bands  (e.g.  shadowing  effects)  while 
the  latter  factor  may  or  may  not  be  correlated,  depending  upon  the  material 
type  of  the  particular  objects  being  viewed.  It  was  decided  to  transform 
these  two  bands  into  two  derived  spectral  bands  which  would,  as  nearly 
as  possible,  isolate  illumination  effects  (a  clutter  factor)  from 
reflectance  effects  (an  information  factor). 


Figure  G-l  is  a  scatterplot  of  the  near- inf rared  data  values  of 
1548  pixels  randomly  selected  from  the  scene.  The  values  plotted  in¬ 
dicate  the  number  of  pixels  having  the  particular  data  value  represented 
(a  "+"  equals  one,  an  "x"  ten  or  more).  The  data  form  a  fan-shaped 
pattern  whose  apparent  vertex,  or  origin,  lies  just  off  the  scale  at 
position  (x,y  =  45,  26).  Physically,  this  origin  corresponds  to  the  zero 
radiance  level  in  both  channels.  The  radial  distance  of  any  data  point 
from  this  origin  is  descriptive  of  its  average  albedo  (average  reflectance 
times  local  illumination  for  both  spectral  bands).  For  any  given  material, 
the  effect  of  variable  illumination  (e.g.  shadows)  can  be  expected  to 
produce  data  variability  along  such  a  radial.  Conversely,  materials 
whose  spectral  reflectance  ratios  for  the  two  bands  are  different  can 
be  expected  to  lie  along  different  radials  (i.e.  along  different  straight 
lines  passing  through  the  apparent  origin  at  different  angles).  In 
order  to  isolate,  as  much  as  possible,  the  shadow  effects  (clutter)  from 
the  reflectance  effects  (information)  the  data  values  of  the  two  near 
infrared  bands,  B1  and  B2,  were  transformed  to  define  two  new  data  bands 
called  "Albedo"  and  "Reflectance  Ratio",  vis. 

T  2  2*1^^ 

Albedo  =  [(Bl  -  26) z  +  (B2  -  45)  J 
Reflectance  Ratio  =  80*arctan  f(Bl  -  26) /(B2  -  45)1 


A  scatter  plot  of  the  resulting  transformed  near  infrared  bands  is 
shown  in  Figure  G-2.  It  will  be  noted  that  most  of  the  rather  high 
degree  of  correlation  seen  in  Figure  G-l  is  missing  from  Figure  G-2, 
thus  implying  that  we  have  at  least  approximately  achieved  our  objective 
of  isolating  independent  information  effects  by  means  of  the  transfor¬ 
mation  used. 

Figure  G-3  is  a  scatterplot  of  the  thermal  infrared  spectral  band 
data  values  for  the  same  1548  pixels  previously  used.  A  high  degree 
of  correlation  is  seen  between  these  bands,  as  is  to  be  expected  based 
upon  the  dominance  of  target  temperature  in  each.  For  the  limited 
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FIGURE  G-l.  SCATTERPLOT  OF  NEAR  INFRARED  SPECTRAL  BANDS 
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temperature  range  involved,  the  correlation  present  is  essentially 
linear,  with  variations  along  the  major  axis  of  the  distribution  repre¬ 
senting  (substantially)  temperature  effects,  and  variations  orthogonal 
to  this  representing  spectral  emittance  differences  between  the  two 
bands . 

An  argument  can  be  made  in  favor  of  using  a  polar  coordinate 
transformation  (as  for  the  reflective  bands)  in  order  to  resolve  independent 
temperature  and  spectral  emittance  variability  in  the  data  of  Figure  G-3. 
However,  in  the  interest  of  simplicity,  and  in  view  of  the  limited  dynamic 
range  of  the  data,  a  cartesian  transformation  of  the  thermal  bands  was 
used  instead.  This  provided  derived  features  which  we  call  "Temperature" 
and  "Emittance  Difference",  with  the  former  representing  distance 
along  the  best  fit  regression  line,  and  the  latter  being  distance  dis¬ 
placed  orthogonal  to  such  a  line.  Numerically,  we  produced  the  derived 
channels: 

Temperature  -  0.75517*B4  +  0.65553*B5  -  30 
Emittance  Difference  -  -0.65553*B4  +  0.75517*B5  +  90 

A  scatter  plot  of  the  resulting  thermal  features  is  shown  in  Figure 
G-4.  Again,  the  transformation  has  essentially  removed  correlation  between 
bands,  and  has  thus  maximized  the  independent  information  content  of 
each  in  a  manner  which  is  physically  interpretable  to  a  first  approximation. 

G.2  SPECTRAL/SPATIAL  CLUSTERING  ALGORITHM 

In  pursueing  the  application  of  remote  sensing  to  agricultural 
inventories,  ERIM  has  developed  a  spectral/spatial  clutering  algorithm 
called  "BLOB".  This  algorithm  operates  on  digital  multispectral  images 
in  such  a  way  as  to  map  out  spatially  continguous  areas  of  homogeneous 
spectral  properties  (blobs).  In  agricultural  applications,  this  algorithm 
is  especially  useful  because  the  resultant  blob  structure  defined  very 
closely  matches  the  actual  agricultural  field  structure  present  in  the 
image,  thus  allowing  subsequent  processing  operations  to  deal  with  the 
mean  spectral  values  of  all  pixels  within  a  field  rather  than  with  each 


146 


Temperature  (digital  counts) 


135.00  + 


157.40 


139.30 


103.20 


74.300 


47.000 


♦♦ 


♦  £ 


♦  ♦  ♦  ♦  ♦  ♦ 

♦  ♦  ♦  4  ♦  ♦  ♦  ♦ 

♦  3  ♦  ♦;?♦♦£  ♦  ♦  ♦  ♦  ♦  ♦  ♦ 

♦♦♦23  4  323^2^2  ♦  ♦  ♦ 

2  2*  ♦  ♦425432*#^»^»  * 

♦  ♦♦54  ♦34792e5^223  2 ♦  ♦  * 

►  ♦  £^43SX6£X44SS545^^3  ♦ 

2425335X  3X  5533  #4  ♦  ♦ 

£♦  ♦  ♦  ♦  3  2. 5  X  ?  ♦  X  X 7X3333*2 3 2 3  < 

♦♦  233227X4X2599X9742342  ♦♦ 

2  ♦S  5 23 X X X X X X X X X 7 7 3 7 ♦ 23 3  2  £♦ 

♦  ♦  2  ♦♦5595X7XXXXXXX.3X3X25  £♦  3  ♦ 

♦♦247XXXX9XXXXX3432232  ♦££ 


♦♦  *♦ 
♦ 


3*XXXXXX69XX9x3X572* ♦ 
£♦ 3  3X997X3XXXX5  4X>  33 
♦324XX5XXXXXxXX7X93^  ♦ 

♦  :•  3 3 334 X  X  X  X  X  X 9 X  X  '3  3  ♦  ♦  2  ♦ 

♦  ♦♦3557^XXXXXXXX2^#  ♦ 

£  ♦  ♦  x  4  5  \  -  v  v  x  X  X  X!  X  2  £  ♦ 

♦35XXXXXXX9* 
♦266XXXXXS 
;?  2453344 


♦ 


+ 

♦ - 

52.000 


-+ - + - 

92.200 


-  + - 4. - +  - 

73.600 


- 4. - 

95.200 


34.400 

Emlttance  Difference  (digital  counts) 


- 4. 

106.0' 


FIGURE  G-4.  SCATTERPLOT  OF  TRANSFORMED  DATA  CHANNELS 
DERIVED  FROM  THERMAL  INFRARED  SPECTRAL  BANDS 
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pixel  as  an  Individual  entity.  This  process  has  obvious  benefit  in  terms 
of  both  precision  (noise  reduction  in  averaging)  and  data  volume  (one 
spectral  vector  per  field  rather  than  one  per  pixel). 

Generically,  BLOB  computes,  for  a  given  pixel,  a  distance  function, 
t,  with  respect  to  each  candidate  blob,  j,  with  which  that  pixel  might  be 
associated.  This  distance  function  is  given  by 


'■  ■  *5(  V*jf 


(l  -  S)2,  (p  -  P,)2] 


(G-l> 


where  x^  =  data  value  of  the  pixel  in  i  spectral  band 
L  =  line  number  of  pixel  in  image 

P  =  point  number  of  pixel  in  image 

x.  =  mean  data  value  of  jtl1  blob  in  1th  spectral  band 

— 1J  th 

L,  =  mean  line  number  of  j  blob 

_ J  ^ 

P,  =  mean  point  number  of  j  blob 

Ax.j_  =  a  characteristic  difference  in  i1-  band  data  values 

Afc  =  a  characteristic  difference  in  line  or  point  number 

n  =■  the  number  of  spectral  bands  used 

A  =  a  spectral  weight  factor  scaled  to  define  the 
relative  ratio  of  spectral/spatial  importance. 

If  the  smallest  value  of  t  found  for  a  given  pixel  is  too  large 
(r,  >  t  J,  the  given  pixel  starts  a  new  blob.  Otherwise,  the  pixel  is 

j  lUaX 

assigned  to  that  blob  minimizing  t.,  and  the  mean  parameters  (x  ,  F  ,  F  ) 

■J  j  j 

for  that  blob  are  updated  accordingly. 

The  spatial  portion  of  tj  is  derived  only  from  the  maximum  of  the 
line  or  point  displacements  from  a  blob  center.  This  feature  recognizes 
the  expected  orientation  relationships  between  rectangular  agricultural 
field  boundaries  and  scan  line  directions  in  Landsat  Multispectral  Scanner 
(MSS)  images,  and  tends  to  produce  rectangular  blobs. 
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Specific  and  detailed  tailoring  of  the  spectral/spatial  clustering 
was  beyond  the  scope  of  this  investigation.  However,  a  first  order  modifi¬ 
cation  of  parameter  settings  was  implemented  in  order  to  reflect  the 
Investigative  intent  of  the  present  work.  It  was  felt  that  Insight  into 
the  nature  of  multispectral  signature  characteristics  of  tactical  targets 
would  be  achieved  by  addressing  the  following  question: 

•  Is  there  any  coherent  spectral  feature  of  a 

target  which  could  separate  it  from  its  immediate 
background  so  as  to  produce  a  recognizable  spatial 
shape? 

The  parameter  settings  for  the  BLOB  algorithm  were  set  in  order  to 
develop  an  understanding  of  the  answer  to  this  question.  » 

G.2.1  Spectral  Parameter  Settings 

The  BLOB  algorithm  of  Equation  (G-l)  is  a  total  distance  measure 
produced  by  adding  a  spectral  distance  function  to  a  spatial  distance 
function.  If  we  ignore  the  latter  term,  we  have  a  purely  spectral  clustering 
algorithm  that  joins  together  spectrally  similar  things  no  matter  where 
in  the  image  they  occur.  In  such  a  situation,  spectral  similarity  occurs 
if 


*  U*1 ' 


±1 


max 


-  A 


(G-2) 


i=l  \ 


Since  ratios  are  involved,  we  have  some  freedom  is  assigning  values  to 

the  three  parameters  t  ,  A  and  Ax.  (x . ,  x..  and  n  are,  of  course,  set 

max  l  l  lj 

by  the  data  used) .  We  arbitrarily  set 


v  ”  30.9  for  convenience 
max 


since  that  is  the  value  used  in  our  agricultural  applications.  For 

similar  reasons,  we  set  Ax  equal  to  1/20  of  the  nominal  data  range  in 
th  ^ 

the  i  spectral  channel  (the  nominal  data  range  is  that  range  of  digital 
values  between  the  5Z  and  95Z  points  on  a  cumulative  histogram  of  the 
channel).  It  then  remained  to  select  a  reasonable  value  for  A. 
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If  we  rewrite  Eq.  (G-2)  for  a  single  spectral  band  (n  »  1),  insert 


the  values  for  Ax  and  t  and  rearrange,  we  have 

max  ® 


x  -  x 


1 


Data  Range 


1/2 

(»)(") 


(G-3) 


<  0-278 

/A 

If  A  =  1,  tnen  Eq.  (G-3)  indicates  that  a  pixel  will  automatically  be 
rejected  from  the  cluster  if  its  data  value, x,  differs  from  the  cluster 
mean,  x^,  by  more  than  the  fraction  0.278  of  the  data  dynamic  range. 

This  is  a  rather  generous  range  for  calling  things  "similar".  Conversely, 
if  A  is  made  very  large,  then  very  small  differences  of  x  -  x^  cause 
rejection  and  data  noise  can  be  the  dominant  feature  driving  the  clustering. 
After  some  investigation,  a  value  of  A  ■  4  was  chosen.  This  means  that 
in  the  single  channel  rendition  of  Eq.  (G-3),  a  difference  from  the 
cluster  mean  of  0.278/A^ -  0.139,  or  about  one-seventh  the  dynamic 
range,  will  define  "similarity".  This  is  about  a  factor  of  four  to  six 
above  the  noise  level  of  a  typical  spectral  channel  evaluated. 

G.2.2  Spatial  Parameter  Settings 


The  primary  purpose  of  the  spatial  portion  of  the  clustering  algorithm 
in  this  investigation  was  to  prevent  spectrally  similar  pixels  which  are 
widely  separated  spatially  from  being  grouped  into  the  same  cluster. 
Otherwise,  the  hope  of  seeing  a  vehicle-shaped  cluster  would  have  been 
compromised.  On  the  other  hand,  it  was  also  necessary  to  choose  the  spatial 
parameters  carefully  so  as  not  to  inadvertantly  carve  up  the  image  into 
vehicle-sized  blobs  in  the  absence  of  any  spectrally  defined  boundaries. 

Both  of  the  above  criteria  were  met  by  choosing  the  characteristic 

length,  Afc,  in  Eq.  (G-l)  to  be  equal  to  2.12  pixels  (or  lines).  This 

means  that,  for  the  previously  defined  t  “30.9,  the  maximum  half- 

nidx 

width  of  a  blob  (i.e.  the  maximum  halfwidth  of  the  region  over  which 
spectrally  similar  pixels  will  be  clustered  together)  is  equal  to 


* 

i 
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max  half  width  *  Ai/r 

max 

=  2.12/30.9 
■  11.8  pixels  (lines) 

At  the  data  scale  of  1.6  feet  per  pixel  (line)  this  means  that  no  cluster 

can  be  bigger  than  about  38  feet  in  length  or  width.  Thus,  if  a  vehicle 

has  spectral  features  that  are  different  from  any  background  within  38 

feet  of  it,  it  will  be  isolated  out  as  a  separate  spatial  entity.  Conversely, 

over  the  spatial  extent  of  a  vehicle  (half-length  »  12  feet  *  7.5  pixels) 

the  maximum  effect  of  the  arbitrary  spatial  cut-off  characteristic  is 

no  more  than  f  .7 *A-)  ,  or  12.5  out  of  the  total  allowable  t  =30.9. 

This  leaves  30.9  -  12.5  =  18.4  units  of  our  distance  measure  for  spectral 

effects.  Thus,  over  the  area  of  even  a  large  vehicle,  the  relative 

weighting  of  spectral  information  changes  little,  from  a  value  of  unity 

1/2 

at  the  center,  to  a  value  of  no  less  than  (18.4/30.9)  =  0.77  at  the 

extreme  edge  of  the  vehicle. 

In  summary,  the  effective  result  of  the  spatial  portion  of  the  algorithm 
is  to  create  a  local  neighborhood  of  less  than  40  feet  diameter  in  the 
vicinity  of  a  target.  If,  within  this  neighborhood,  the  target  has  some 
internally  consistent  spectral  signature  which  differs  from  that  of  the 
local  backgrounds  present,  the  BLOB  spatial/spectral  clustering  algorithm 
will  isolate  a  target  shaped  entity  (blob)  which  may  subsequently  be 
recognized  by  its  spatial  features. 

G.3  PROCESSING  AND  RESULTS 

The  BLOB  spectral/spatial  clustering  algorithm  was  used  with  the 
derived  spectral  features  defined  in  Section  G.l  in  order  to  assess  the 
extent  to  which  there  exists  a  coherent  spectral  signature  for  a  target 
vehicle  which  allows  it  to  be  isolated  from  its  immediate  background 
to  produce  a  recognizable  spatial  shape.  The  derived  spectral  feature 
channels  for  the  1330  hours  Camp  A.P.  Hill  data  were  utilized  in  various 
combinations  as  shown  in  the  test  cases  defined  in  Table  G-l. 
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TABLE  G-l.  SPECTRAL/SPATIAL  CLUSTERING  CHANNEL  COMBINATIONS 


Test  Case 


* 

Derived  Spectral  Feature  Channels 


1 


Albedo,  Reflectance  Ratio,  Temperature, 
Emittance  Difference 


2  Reflectance  Ratio,  Temperature, 

Emittance  Difference 

3  Reflectance  Ratio,  Emittance  Difference 

4  Temperature,  Emittance  Difference 

5  8.4  to  13.7  vim  (original  spectral  channel) 


* 

See  Section  G.l. 


* 


52 


Test  Case  1  used  all  of  the  derived  spectral  features  available 
and  represented  the  maximal  use  of  all  independent  spectral  information 
present  in  the  original  multlspectral  data.  Test  Case  2  dispensed  with 
the  albedo  feature,  thus  providing  for  a  first  order  elimination  of  the 
clutter  effects  caused  by  shadowing.  Test  Case  3  further  eliminated 
temperature  and  thus  used  only  those  features  associated  most  directly 
with  the  intrinsic  optical  properties  of  the  scene  materials  present. 

Test  Case  4  included  only  the  independent  thermal  information  features 
provided  by  the  two  thermal  infrared  spectral  bands  originally  collected. 
Finally,  Test  Case  5  utilized  only  a  single  thermal  infrared  spectral 
band,  one  of  the  original  spectral  bands  collected. 

Spectral/spatial  clustering  was  performed  on  the  five  test  cases 
using  the  1330  hours  Camp  A.P.  Hill  data  base.  Figure  G-5  shows  a  small 
portion  of  an  output  map  for  one  of  the  test  cases.  Each  blob  isolated 
is  keyed  by  a  different  symbol  combination.  Since  there  are  literally 
thousands  of  blobs  produced  in  each  test  case,  a  map-like  display  is 
somewhat  of  a  problem  in  terms  of  coding.  Me  have  not  very  effectively 
solved  that  problem,  and  any  display  more  comprehensive  than  that  of 
Figure  G-5  rapidly  becomes  unintelligable  in  terms  of  presentation  in  a 
report  such  as  this. 

Analysis  of  the  processed  results  was  designed  to  proceed  in  four 
serial  steps: 


1) 

Determination  of  whether  a  target  was  isolated  as  a 

vehicle  shaped  blob; 

2) 

If  1)  was  successful,  definition  of  a  spatial 

filtering  algorithm/approach  to  cueing  those  blobs; 

1 

3) 

If  1)  and  2)  were  successful,  assessment  of  the  number 
of  non-target  blobs  also  cued  (false  alarms); 

- 

and  4) 

If  false  alarms  from  3)  were  substantial,  investigation 

to  determine  whether  additional  false  alarm  rejection 

could  be  accomplished  by  some  spectral  selection  criterion 

following  spatial  detection. 
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FIGURE  G-5.  EXAMPLE  DISPLAY  OF  SPECTRAL/SPATIAL  CLUSTERS, 
TEST  CASE  1  INDIVIDUAL  BLOB  PIXELS  ARE  SYMBOL-KEYED 


Unfortunately,  the  analysis  results  from  Step  1  quickly  highlighted  the 
fact  that  serious  problems  existed  in  the  basic  concept  that  there  might 
be  something  coherent  about  the  spectral  signatures  within  an  Individual 
target.  In  general,  it  was  found  that  most  of  the  targets  present  were 
actually  represented  by  several  very  small  blobs.  Apparently,  the  within- 
target  variance  of  all  feature  combinations  tried  was  just  too  great  in 
general  to  create  a  single  major  spatial  feature  possessing,  by  itself, 
any  useful  shape  information. 

Table  G-2  provides  a  summary  of  the  analysis  results  partitioned  by 
test  case  and  by  target.  Some  exceptions  to  the  general  negative  result 
can  be  seen.  In  particular,  the  3/4  ton  truck  and  the  5  ton  truck  generally 
were  isolated  out  of  their  backgrounds  by  virtue  of  creating  vehicle-sized 
rectangular  blobs.  On  the  other  hand,  one  of  the  APC's  and  the  2-1/2  ton 
truck  were  never  isolated  from  their  immediate  backgrounds,  not  even  as  a 
grouping  of  small  blobs. 

The  results  obtained  do  identify  some  relative  utility  to  the  derived 
spectral  features  utilized.  The  best  results,  both  overall  and  for  each  tar¬ 
get,  were  obtained  in  Test  Case  2.  This  supports  the  original  hypothesis  by 
the  shadowing  effects  represented  by  the  albedo  feature  which  was  eliminated 
in  Test  Case  2.  On  the  other  hand,  the  worst  results  were  found  in  Test 
Case  3.  Apparently,  there  is  not  a  large  amount  of  discriminative  power 
in  the  surface  optical  properties  alone,  at  least  as  represented  by  the 
specific  feature  transformations  used  on  the  data  set  processed. 

The  implications  of  the  results  achieved  are  significant.  We  have 
taken  one  approach  to  transforming  the  raw  spectral  data  into  derived 
features  which  should  be  more  stable,  and  certainly  are  more  intuitively 
understandable.  We  have  asked  whether  there  is  any  uniformity  of  those 
features  on  a  tactical  target,  under  conditions  where  the  term  "uniformity" 
represents  a  rather  substantial  allowance  for  variance  with  respect  to 
the  dynamic  range  of  the  data.  And  we  have  found  in  general,  and  in 
every  instance  for  armored  vehicles,  that  the  answer  to  that  question, 
at  least  for  this  limited  data  set, is  negative. 
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TABLE  G-2 .  SPECTRAL/SPATIAL  CLUSTERING  RESULTS 


Target 

Location 

Spectral/Spatial 

1 

Cases 

K 

(Line/Point) 

r 

2 

3 

4 

~5 

APC 

187/443 

P 

P 

P 

P 

P 

2  1/2  T 

228/481 

- 

- 

- 

- 

- 

M-60 

303/414 

P 

P 

- 

- 

- 

M-60 

343/511 

P 

P 

- 

P 

P 

3/4  T 

388/513 

P 

/ 

- 

/ 

/ 

APC 

406/422 

- 

? 

- 

- 

- 

5  T 

435/463 

/ 

✓ 

/ 

/ 

✓ 

APC 

445/430 

- 

- 

- 

- 

- 

Key:  / 

=  vehicle  shaped  blob 

isolated 

P 

=  several  small  blobs 

isolated 

which  add 

up  to 

a  / 

- 

=  target  not 

isolated 

from  background  in 

any  manner 

* 

See  Table  G-l. 


156 


REFERENCES 


[1]  Stewart,  Steve,  MS  Cueing  Techniques  Program,  Mission  Report  for 
Night  Vision  Laboratory  Program,  Camp  A.P.  Hill  Virginia  Test 
Site,  March  28,  29,  30,  1978,  Environmental  Research  Institute  of 
Michigan,  Ann  Arbor,  Mi.  Report  No.  131500-8-P,  May  1978. 

[2]  Spectral  and  Polarization  Characteristics  of  Selected  Targets  and 
Backgrounds:  Instrumentation  and  Measured  Results  93.3-14.0  pm), 

D.  Faulkner,  R.  Horvath,  J.P.  Ulrich,  E.  Work,  Willow  Run  Labora¬ 
tories,  The  University  of  Michigan,  Technical  Report  No.  AFAL-TR- 
71-199,  August  1971. 

[3]  Handbook  of  Military  Infrared  Technology,  William  L.  Wolfe,  Editor, 
The  University  of  Michigan,  1965,  Sponsored  by  Office  of  Naval 
Research,  Department  of  the  Navy,  Washington,  D.C. 

[4]  M.  Bair,  et  al.,  "Active  Passive  Cueing  Techniques"  (U) ,  ERIM 
Report  No.  112900-69-F,  Air  Force  Avionics  Laboratory,  Contract 
No.  F33615-75-C-1155,  (In  Publication  June  1979)  (SECRET). 

[5]  S.  Sternberg,  et  al.,  "3D  Target  Classifier"  (U),  ERIM  Report  No. 
126200-36-F,  Air  Force  Avionics  Laboratory,  Contract  No.  F33615-77- 
C-1010  (In  Publication  June  1979)  (CONFIDENTIAL) 

[6]  P.G.  Hasell,  Jr.,  et  al.,  "Michigan  Experimental  Multispectral 
Mapping  System,  A  Description  of  the  M7  Airborne  Sensor  and  Its 
Performance,"  Report  No.  190900-10-T,  Environmental  Research 
Institute  of  Michigan,  Ann  Arbor,  Michigan,  January  1974. 

[7]  Hasell,  P.  G.,  and  L.  M.  Larsen,  "Calibration  of  an  Airborne  Multi¬ 
spectral  Optical  Sensor,  Report  6400-137-T,  Willow  Run  Laboratories, 
Institute  of  Science  and  Technology,  The  University  of  Michigan, 

Ann  Arbor,  Michigan,  September  1968. 


DISTRIBUTION  LIST 


W26P70  10 

Night  Vision  Laboratory 
USAECOM,  DRSEL-NV-VI 
Fort  Belvoir,  VA  22060 

Defense  Documentation  Center  14 

ATTN:  DDC-TCA 

Cameron  Station  (Bldg  5) 

Alexandria,  VA  22314 

Director  1 

National  Security  Agency 
ATTN:  TDL 

Fort  George  G.  Meade,  MD  20755 

Advanced  Research  Projects  Agency  1 

1400  Wilson  Blvd 
Arlington,  VA  22209 

Conmander,  Naval  Ship  Systems  Cmd  1 

Technical  Library,  Rm  3  S08 
National  Center  No.  3 
Washington,  DC  20360 

Department  of  the  Navy  1 

Chief,  Office  of  Naval  Research 
Washington,  DC  20350 

Director  1 

Naval  Research  Laboratory 
Code  2627 

Washington,  DC  20390 

Conmander  1 

Naval  Electronics  Laboratory  Ctr 
ATTN:  Library 
San  Diego,  CA  92152 

Conmander  1 

U.S.  Naval  Ordnance  Lab 
ATTN:  Technical  Library 
White  Oak,  Silver  Spring,  MD  20910 

Conmander  2 

Naval  Electronic  Systems  Cmd  Hq 
ATTN:  Code  0563 
Washington,  DC  20360 

Conmandant,  Marine  Corps  1 

Hq. ,  U.S.  Marine  Corps 
ATTN:  Code  A04C 
Washington,  DC  20380 

Communications  Electronics  Div.  1 

Development  Center 
Marine  Corps  Dev  &  Educ  Cmd 
Quantico,  VA  22134 


Marine  Corps  Ln  Officer  1 

Bay  2C214A,  Hex  Bldg 
ATTN:  John  W.  Everett 
China  Lake,  CA  93555 

Naval  Weapons  Center  1 

Code  5525 

ATTN:  Henry  Blazak 
China  Lake,  CA  93555 

Rome  Air  Development  Center  1 

ATTN:  Documents  Lib  (TDLD) 

Griff iss  AFB ,  NY  13440 

Hq.  ESD  (TRI)  1 

L.G.  Hanscom  Field 
Bedford,  MA  01730 

Air  Force  Avionics  Laboratory  2 

ATTN:  AFAL/DOT,  STINFO 
Wright-Patterson  AFB,  OH  45433 

Recon  Central/RSA  1 

Air  Force  Avionics  Laboratory 
Wright-Patterson  AFB,  OH  45433 

Air  Force  Avionics  Laboratory  1 

ATTN:  AVRO 

Wright-Patterson  AFB,  OH  45433 

Armament  Dev  &  Test  Center  1 

ATTN:  SSLT 

Eglin  Air  Force  Base,  FL  32542 

Hq.  Air  Force  Systems  Cmd  1 

ATTN:  DLTE 
Andrews  Air  Force  Base 
Washington,  DC  20331 

Tactical  Air  Reconnaissance  Ctr  1 
Army  Liaison  Office 
Shaw  AFB.  SC  29152 

HQDA  (DAMI-ZA)  2 

Washington,  DC  20310 

Commander  1 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATTS-X 
Fort  Monroe,  VA  23651 

HQDA  (DACE-ED)  1 

Washington,  DC  20314 

Conmander  1 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCD-SI 
Fort  Monroe,  VA  23651 


159 


HQDA  (DARD-DDC)  1 

Washington,  DC  20310 

Commander 

U.S.  Army  Training  &  Doctrine  Cmd  1 

ATTN:  ATCD-C 

Fort  Monroe,  VA  23651 

HQDA  (DARD-ARP/Dr.  R.  B.  Watson)  1 

Washington,  DC  20315 

Commanding  General  1 

U.S.  Army  Materiel  Command 
ATTN:  AMCMA-EE 
Washington,  DC  20315 

Commanding  General  1 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-FW 

Washington,  DC  20315 

Commander  1 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCD-F 
Fort  Monroe,  VA  23651 

Commanding  General  1 

U.S.  Army  Materiel  Command 
ATTN:  AMCRD-H 
Washington,  DC  20315 

Commanding  General  1 

U.S.  Army  Missile  Command 
ATTN:  AMCMI-RR,  Dr.  J.P.  Hallowes 
Redstone  Arsenal,  AL  35809 

Conmanding  General  2 

U.S.  Army  Missile  Command 
Redstone  Scientific  Info  Center 
ATTN:  Chief,  Document  Sec 
Redstone  Arsenal,  AL  35809 

Commander  1 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCE 
Fort  Monroe,  VA  23651 

Commanding  General  1 

U.S.  Army  Weapons  Command 
ATTN:  AMSWE-REF 
Rock  Island,  IL  61201 

Commanding  General  1 

U.S.  Army  Weapons  Command 
ATTN:  AMSWE-RER-L 
Rock  Island,  IL  61201 

Commander  1 

U.S.  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCAIC 

Fort  Leavenworth,  KS  66027 


Commander  1 

U.S.  Army  Training  &  Doctrine  Cmd 
ATTN:  ATCD-CI 
Fort  Monroe.  VA  23651 

Comnander  1 

U.S.  Army  Logistics  Center 
ATTN:  ATCL-MC 
Fort  Lee,  VA  23801 

Commandant  1 

U.S.  Army  Ordnance  School 
ATTN:  ATSOR-CTD 

Aberdeen  Proving  Ground,  MD  21005 

Commander  1 

U.S.  Army  Intelligence  School 
ATTN:  ATSIT-CTD 
Fort  Huachuca,  AZ  85613 

Commandant  1 

U.S.  Army  Field  Artillery  School 
ATTN:  ATSFA-CTD 
Fort  Sill,  OK  73503 

Commandant  1 

U.S.  Army  Air  Defense  School 
ATTN:  ATSAD-CTD 
Fort  Bliss,  TX  79916 

Commandant  1 

U.S.  Army  Engineer  School 
ATTN:  ATSEN-CTD 
Fort  Belvoir,  VA  22060 

Commanding  Officer  1 

Harry  Diamond  Laboratories 
ATTN:  Library 
Washington,  DC  20438 

CO,  USA  Foreign  Science  &  Tech  Ctr  2 

ATTN:  AMXST-ISI 
220  Seventh  St,  NE 
Charlottesville,  VA  22901 

CO,  USA  Foreign  Science  Div  1 

ATTN:  AMXST  CE  Division 
220  Seventh  St,  NE 
Charlottesville,  VA  22901 

CO.  US  Army  Picatinny  Arsenal  1 

ATTN:  SMUPA-VC5  (Mr.  P.  Kisatsky) 

Bldg  350 

Dover,  NJ  07801 

CO,  US  Army  Picatinny  Arsenal  2 

ATTN:  SMUPA-RT-S 
Bldg  59 

Dover,  NJ  07801 

ConmandinR  Officer  1 

Frankford  Arsenal 
ATTN:  Library,  H1300.  Bl.51-2 
Philadelphia,  PA  19137 


160 


2 


Conroanding  Officer  1 

Frankford  Arsenal 
ATTN:  W1000-65-1  (Mr.  Helfrlch) 
Philadelphia,  PA  19137 

Commanding  Officer  1 

Frankford  Arsenal 
ATTN:  SMUFA-W- 1000  (Mr.  Kerensky) 
Philadelphia,  PA  19137 

Commanding  General  1 

White  Sands  Missile  Range 

ATTN:  STEWS- ID- S  Hq 

White  Sands  Missile  Range,  NM  88002 

Commanding  Officer  1 

Edgewood  Arsenal 
ATTN-  SMUEA-TSTI-TL 
Edgewood  Arsenal,  MD  21010 

Conmanding  Officer  1 

U.S.  Army  Mtls  &  Mech  Rsch  Center 
ATTN:  AMXMR-ATL ,  Tech  Lib  Br 
Watertown,  MA  02172 

President  1 

U.S.  Army  Artillery  Board 
Fort  Sill,  OK  73503 

Commanding  Officer  1 

Aberdeen  Proving  Ground 
ATTN:  STEAP-TL 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

Aberdeen  Proving  Ground 
ATTN:  STEAP-MT-TF,  Bldg  436 
Aberdeen  Proving  Ground,  MD  21005 

Director,  Ballistic  Rsch  Labs  1 

U.S.  Army  Aberdeen  R&D  Center 
ATTN:  AMXRD-BTL  (Mr.  F.  J.  Allen) 
Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

U.S.  Army  Ballistic  Rsch  Labs 
ATTN:  AMXRD-BSP  (Dr.  D.  Reuyl) 
Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

U.S.  Army  Land  Warfare  Lab 
ATTN:  CRDLWL-1A 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

U.S.  Army  Limited  Warfare  Lab 
ATTN:  CRDLWL-7C 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

U.S.  Army  Land  Warfare  Lab 
ATTN:  Mr.  David  Samuels 
Aberdeen  Proving  Ground,  MD  21005 


Commanding  Officer 

U.S.  Army  Electronic  Proving  Grd 

ATTN:  STEEP-MT 

Fort  Huachuca.  AZ  85613 


Chief,  USAECOM  Field  Engr  Ofc  1 

U.S.  Army  Electronic  Proving  Grd 
ATTN:  STEEP-LN-F  (Mr.  H.A.  Ide) 

Fort  Huachuca,  AZ  85613 

Commanding  Officer  1 

USASA  Test  &  Evaluation  Center 
Fort  Huachuca,  AZ  85613 

U.S.  Army  Research  Ofc-Durham  1 

ATTN:  CRDARD-IP 
Box  CM,  Duke  Station 
Durham,  NC  27706 

U.S.  Army  Research  Ofc-Durham  1 

ATTN:  Dr.  Robert  J.  Lontz 
Box  CM,  Duke  Station 
Durham,  NC  27706 

USA  Security  Agcy  Combat  Dev  Act.  1 
ATTN:  IACDA-EW 

Arlington  Hall  Station,  Bldg  420 
Arlington,  VA  22212 

Commandant  1 

U.S.  Army  Air  Defense  School 
ATTN:  C&S  Dep  MSL  Sci  Div 
Fort  Bliss,  TX  79916 


Commander  * 

U.S.  Army  Combined  Arms  Combat 
Developments  Activity 
ATTN:  ATCACC 

Fort  Leavenworth,  KS  66027 

Commanding  Officer  1 

U.S.  Army  Dugway  Proving 
Ground  Lib 

ATTN:  STEDP-TL,  Tech  Lib 
Dugway,  UT  84022 

Commanding  Officer  1 

U.S.  Army  Cold  Regions  R&E  Lab 
ATTN:  CRREL-RP  (Dr.  Yin  Chao  Yen) 
Hanover,  NH  03755 

Commanding  Officer  1 

Yuma  Proving  Ground 

ATTN:  STEYP-AD  (Tech  Lib) 

Yuma,  AZ  85364 

Director  1 

U.S.  Army  Advanced  Mtl  Concepts  Agcy 
ATTN:  AMXAM 
Washington,  DC  20315 

Commanding  Officer  1 

U.S.  Army  Materiel  Comnand 
ATTN:  AMCRD-R  (H.  Cohen) 

Washington ,  DC  20315 


161 


1 


1 


Commandant 

U.S.  Army  Infantry  School 
ATTN:  ATSIN-CTD 
Fort  Benning,  GA  31905 

Commander  1 

U.S.  Army  Logistics  Center 
ATTN:  ATCL-X 

Aberdeen  Proving  Ground,  MD  21005 

Commandant  1 

U.S.  Army  Armor  School 
ATTN:  aTSAR-CTD 
Fort  Knox,  KY  40121 

Commandant  1 

U.S.  Army  Aviation  School 
ATTN:  ATSAV-CTD 
Fort  Rucker,  AL  36360 

Commandant  1 

U.S.  Army  Institute  for  Mil  Assistance 
ATTN:  ATSMA-CTD 
Fort  Bragg,  NC  28307 

Comnandant  2 

U.S.  Army  Field  Artillery  School 
ATTN:  Target  Acquisition  Dept 
Fort  Sill,  OK  73503 

Commanding  General  1 

U.S.  Army  Missile  Conmand 
ATTN:  AMSMI-RFG  (Mr.  N.  Bell) 

Redstone  Arsenal,  AL  35809 

Commanding  Officer  1 

Harry  Diamond  Laboratories 
ATTN.  AMXDO-RCB  (Mr.  J.  Nemarich) 
Washington ,  DC  20438 

Commanding  Officer  1 

U.S.  Army  Systems  Analysis  Agcy 
ATTN:  AMSRD-AMB  (Mr.  A.  Reid) 

Aberdeen  Proving  Ground,  MD  21005 

Commanding  Officer  1 

Frankford  Arsenal 
Joint  Laser  Safety  Team 
ATTN:  SMUFA-W1000 
Philadelphia,  PA  19137 

Commanding  General  1 

Project  MASSTER 
Fort  Hood,  TX  76544 

U.S.  Army  Liaison  Office  1 

MIT-Lincoln  Lab,  Rm  A-210 
P.0.  Box  73 
Lexington,  MA  02173 

Chief  1 

Missile  Electronic  Warfare  Tech  Area 
EW  Lab,  USA  Electronics  Conmand 
White  Sands  Missile  Range,  NM  88002 


CG,  U.S.  Army  Electronics  Cmd 
ATTN:  AMSEL-PP/P-IED 

(Mr.  C.  Mogavero) 

225  South  18th  St 
Philadelphia,  PA  19103 

Chief,  Intelligence  Materiel  1 

Development  Office 
Electronic  Warfare  Lab, 

USAECOM 

Fort  Holabird,  MD  21219 

Commanding  General 

U.S.  Army  Electronics  Command 

Fort  Monmouth,  NJ  07703 


AMSEL-WL-D  1 

AMSEL-NL-D  1 

AMSEL-TL-D  1 

AMSEL-VL-D  1 

AMSEL-BL-D  1 

AMSEL-CB  1 

AMSEL-SI-CB  1 

AMSEL-MA-MP  1 

AMSEL-MS-TI  2 

AMSEL-GG-TD  1 

CSCS-FM  1 

AMSEL-RD  1 

TRADOC-LNO  1 

Director 

Night  Vision  Laboratory 

Fort  Be 1 voir ,  VA  22060 

AMSEL-NV-D  1 

AMSEL-NV-VI  1 

AMSEL-NV-PA/RA  1 

Commanding  General  1 


U.S.  Army  Strategic  Conmunications 
Command 

ATTN:  SCC-FM-D 

Fort  Huachuca,  AZ  85613 

Sylvania  Electronic  Systems  1 

Western  Division 
ATTN:  Tech  Reports  Library 
P.  0.  Box  205 
Mountain  View,  CA  94040 

NASA  Scientific  &  Tech  Info  Facility  2 
ATTN:  Acquisitions  Br  (S-AK/DL) 

P.  0.  Box  33 

College  Park,  MD  20740 

Advisory  Group  on  Electron  Devices  2 
ATTN:  Secy,  SP  Gr  on  Optical  Masers 
201  Varick  St 
New  York,  NY  10014 

Electronic  Properties  Info  Ctr  1 

Hughes  Aircraft  Company 
Centinela  and  Teale  Streets 
Culver  City,  CA  90230 

Remote  Area  Conflict  Info  Ctr  1 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  OH  43201 


